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ARTIST'S CONCEPTION OF LUNAR WORM 


FEASIBILITY STUDY FOR LUNAR WORM 


PLANETARY ROVING VEHICLE CONCEPT 

By. F. A. Dobson and D. G. Fulton 
Aeronutronic, Division of Philco Corporation 

SUMMARY 


A feasibility study has been carried out of various concepts of using a 
bellows to provide cross country mobility, particularly on the moon. Effort 
was directed primarily toward analytical treatment of bellows mobility and 
of bellows structure including both corrugated and membrane bellows, and 
applying the results to several brief design studies of typical missions. All 
concepts studied can have all moving parts internally located in a protected, 
pressurized environment and can be considerably foreshortened during 
delivery. 

The concept wherein the bottom surface of the vehicle moves in the form 
of a traveling wave similar to the motion used by snakes and centipedes 
appears to be potentially the most versatile.* As compared with current cross 
country vehicles, it is capable of moderate to high speeds over rough terrain 
with a smoother ride and less vehicle bounce, good maneuverability, good 
propulsive efficiency, better ability to climb obstacles and to bridge crevices. 
It can be designed with a very low footprint pressure for use over a wide 
variety of surfaces including even fluid-like soils through which it is capable 
of swimming. 

The double-acting bellows concept, which is a variation of rib-walking, 
where adjacent ribs contacting the surface are 180° out of phase, appears 
quite suitable for unmanned vehicles with regard to simplicity, ease of 
deployment, weight and power requirements but the ride tends to be rough 
and ability to clear obstacles is somewhat limited. A version of this concept 
can be designed which is within the present state-of-the-art of metal bellows 
technology. 

The extension-contraction bellows concept, wherein the vehicle alternately 
shortens and lengthens and propulsion is derived from asymmetry built into 
the lower surface, is the simplest concept considered but also requires the 
most power because of friction in sliding over the terrain. It is restricted in 
footprint loading and maneuverability but is uniquely suited for incorporation 
of a novel solar-mechanical energy transfer system which holds some 
promise as a lightweight, inexpensive means for supplying mechanical energy 
for propulsion. 


The mechanization of the traveling wave motion is the subject of a prior, 
independent, Aeronutronic invention disclosure dated February 5, 1965, 
entitled "Peristaltic Track. " 


The study also includes fundamental analytical treatment of a stepped 
surface planing over soil, of wave motions which might be used by a variety 
of vehicles and of the structure of membrane bellows. A lunar soil model 
suitable for short stroke and cyclic motion studies was developed. 
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INTRODUCTION 


An important part of the overall program for lunar exploration is trans- 
portation over the surface of the moon. The characteristics of this surface 
are still largely unknown and are the subject of widely differing predictions. 

It is agreed, however, that in certain respects (such as reduced gravity, 
extreme vacuum, severe temperature extremes, and various hazards such 
as radiation and meteorites), the environment will differ greatly from that on 
earth. For these reasons, lunar vehicles undoubtedly will have to possess 
new and unusual features. These will be in addition to the basic off- road 
capability which is necessary in traversing unprepared surfaces. A variety 
of concepts for providing this mobility have been considered, some of which 
depart substantially from conventional wheeled and tracked vehicles such as 
walkers, walking-beam wheels, elastic wheels, and jumpers (pogo type). 

The purpose of this study is to investigate the possibility of satisfying 
the special requirements for lunar vehicles by use of new concepts based on 
sealed expandable bellows. A bellows offers significant advantages including 
no external moving parts in the harsh lunar environment, low footprint pres- 
sure desirable for surfaces of unknown texture or bearing strength, ability 
to bridge crevices, convenience for packaging, and possibility of increasing 
effective wall thickness for emergency protection from radiation during solar 
flares. Because of the similarity between an elongated bellows and various 
types of worms or caterpillars, it is logical to investigate the use of worm- 
like motions of the bellows for propulsion. For this reason, this concept is 
named the "Lunar Worm. " The major task of the study is to determine the 
feasibility of a bellows for providing vehicle mobility on the lunar surface 
with secondary tasks to explore in a preliminary manner its use for several 
possible mission applications, to provide a gross comparison between the 
performance of bellows vehicles and other type vehicles and to point out the 
direction of further effort. 

Mobility analyses were made for extension-contraction type bellows 
vehicles, for rib-walking type bellows vehicles, and for traveling wave type 
bellows vehicles. Structural analyses were made for preformed bellows and 
inflatable membrane bellows including the development of an analytical 
technique for the latter. Design studies were made for movable lunar shelters, 
a small unmanned vehicle having limited range and a manned roving vehicle 
having extended range. Overall evaluation was made from considerations of 
mobility, structural limitations, and propulsive energy. 

Information on the present state-of-the-art in convoluted metal bellows 
design was furnished by Solar Division of International Harvester Corporation, 
under subcontract to Aeronutronic. The design study of the Small Unmanned 
Vehicle is based upon designs proposed by Solar. 
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SYMBOLS 


A 


a 

c 

a 

e 

B 


b 


C 

c 

D 

D 

e 

D /W 
e 

d 

E 


F 

f 

G 

g 

g e 

H 

h 

I 

J 


Amplitude of a wave; length of the cranks of a traveling-wave 
device 

Contracted length of a vehicle segment 
Extended length of a vehicle segment 

Location of the point on a wave surface that is fixed with respect 
to the ground 

Width of vehicle; width of track; for bellows design, the effec- 
tive vehicle width or the distance between points where bending 
moment changes sign 

Constant of proportionality for ELMS model 
Coefficient of soil cohesion 
Drag force; external drag; diameter 
Equivalent drag 

Equivalent drag-to-weight ratio 
Bellows depth 

Modulus of elasticity; location of the intersection of a wave 
surface with the ground surface 

6/L vehicle strain; with an exponent - the Naperian base 

Allowable stress; allowable fatigue stress 

Frequency; internal stress 

Modulus of rigidity 

Acceleration due to local gravity 

Acceleration due to Earth 1 s gravity 

Tractive force 

Vehicle height 

Moment of inertia 
2 

(j! y )/m dimensionless planing parameter 

o ir 
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Soil displacement from its initial location 

Modulus of soil deformation due to cohesion 

Modulus of soil deformation due to friction 

Soil parameter reflecting degree of soil compactness 

Soil parameter reflecting the fundamental character of the soil 
shear curve 

Effective coefficient of friction of the aft-moving portion of an 
extension- contraction vehicle 

Effective coefficient of friction of the forward— moving portion of 
an extension-contraction vehicle 

Closed (stored) length of a bellows vehicle 

Vehicle length; of an extension-contraction vehicle - the minimum 
length during an operating cycle 

Length of stroke of an extension- contraction vehicle 

Effective length in the direction of motion of a planing surface 

Mass; total vehicle mass 

Bending moment 

Empty vehicle mass 

Mass of useful load 

Loading on a planing surface per unit of width normal to the 
direction of motion 

Number of bellows convolutions in a vehicle 
Number of segments of a segmented vehicle 

Number of segments of a vehicle which are extended; exponent, 
dimensionless factor reflecting stratification of the soil 

Force; force exerted on the end of the vehicle 

Bellows spring force 

Frictional force 

Power required to move a wave surface on a moving vehicle; 
power required to propel a vehicle 



Point on an extension- contraction vehicle which is instantane- 
ously at rest with respect to the ground; vehicle footprint 
pressure 

Footprint pressure of empty vehicle 
Vehicle planform loading 

Maximum allowable planform loading for a vehicle to plane 
according to the planing theory 

Internal pressure in a vehicle 

Running load in a member per unit of length (width) 

Ratio of vehicle length per foot (planing surface) to the effective 
length of the foot in the direction of motion 

Spring rate 

Radius (of vehicle); hoop radius; horizontal resistance to motion 
along the shear plane 

Area; footprint area 

Planform area 

Distance traveled 

Period, time per cycle 

Time; bellows material thickness 

Effective thickness 

Volume; vehicle enclosed volume 

Vehicle (average) forward velocity 

Tangential velocity of a wave surface with respect to the vehicle 

Wave velocity relative to the ground 

Weight in the local gravity; total vehicle weight 

Weight of one segment of a segmented vehicle 

Location of a point on a wave surface 

Distance between the node of the transverse wave and the 
corresponding node of the longitudinal wave 



*X 

X, y, z 

y 


a 

$ 

y 

6 

A 


V 

9 

\ 


P 

*a 


V 

i 

ir 

P 

a 

T 


Longitude coordinate of a wave 
Generalized orthogonal coordinate 

Distance from the neutral axis to the extreme fiber (in calculat- 
ing bending stress) 

Depth of penetration of the soil surface 

Angle of attack of a planing surface with respect to the vehicle 

Slope of a wave surface 

Specific mass of structural material 

Specific mass of the soil 

Bellows extension per disc 

Incremental value 

Efficiency 

Slope angle 

Wavelength 

Friction coefficient 

Soil friction coefficient; tan <p 

Coefficient of friction of the aft-moving portion of an extension- 
contraction vehicle 

Coefficient of friction of the forward-moving portion of an 
extension-contraction vehicle 

Poisson ! s ratio 

Particle displacement 

3 . 14 

Local radius of curvature 
Normal soil pressure 
Soil shear stress 
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I 

1 

I, 

<p Angle of friction between soil grains 

0 Phase angle between adjacent cranks of a traveling wave vehicle I 

00 Angular velocity j 
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MOBILITY ANALYSIS 


In this section, generalized expressions are developed for velocity and 
power requirements of various configurations as a function of soil charac- 
teristics. Although a number of simplifying assumptions are made, the 
results are considered accurate enough for use in carrying out gross com- 
parisons between different systems. 

The statement of work for this study required the use of the NASA Engi- 
neering Lunar Model Surface (ELMS). It subsequently developed that ELMS 
had been developed for wheeled or tracked vehicle analyses but was not 
adequate for the present study which required analysis involving cyclic motion 
of the vehicle contacting the surface. In addition, certain of the ELMS rela- 
tions were dimensionally inconsistent for analytical studies of this nature. 
Therefore, it became necessary to develop soil models which were both 
appropriate and realistic for the present study. 

To obtain a clear picture of the mobility mechanism, the bottom surface 
of the vehicle which actually contacts the soil is treated as a thrust-producing 
mechanism. All sources of drag, such as bulldozing, sliding friction along 
the sides of the vehicle, climbing angle, and drag due to some object or 
vehicle being towed, are lumped together as ’’external drag. 11 The ratio of 
external drag to total vehicle weight is treated as an equivalent climbing 
slope, which is considered to be the independent variable in most of the charts 
which follow. Since the traction is obtained by friction with the soil, the 
equivalent climbing slope can never be greater than the average vehicle-to- 
soil friction coefficient. The velocities are handled in nondimensional form. 

The power required for propulsion under various conditions is given by 
an ’’equivalent drag-to-weight ratio, ” which is equal to the propulsive power 
divided by the vehicle weight and the average velocity. This ratio is called 
"specific power” by some authors. For a propulsive system with no losses, 
the equivalent drag-to-weight ratio would be equal to the equivalent slope. In 
an actual case, the equivalent drag-to-weight will always be greater than the 
equivalent slope and will become infinite when the vehicle stalls. When the 
equivalent slope becomes negative (that is, on a down grade) the equivalent 
drag-to-weight becomes quite small, and in some cases may become nega- 
tive indicating that power may be recovered from the propulsion system. 

The mobility studies are divided into three major sections, the first of 
which is concerned with the Extension- Contraction motion, the second with 
the Rib-Walking motion and the third with Traveling Wave motion, which is 
a generalization of rib walking. 
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LUNAR SOIL MODEL 


For quantitative comparisons of different mobility systems, it is neces- 
sary to set up a soil model. Since the NASA Engineering Lunar Model Surface 
(ELMS), as mentioned previously, was found to have shortcomings, two other 
models were introduced. In this section the three models will be discussed 
briefly. A more detailed description of the ELMS will be found in Appendix B. 


In an analysis of vehicle mobility, the soil characteristics of greatest 
importance is the friction coefficient between the vehicle traction surface and 
the soil. This friction coefficient, is equal to T/cr, the ratio of shear 
stress to normal stress at the interface between the soil and the vehicle. The 
value of T/ct specified in the ELMS is shown as a function of soil slippage, or 
horizontal displacement in curve III, Figure 1. 



SOIL SLIPPAGE, j (INCHES) 


F05396 U 


FIGURE 1. SOIL FRICTION COEFFICIENT VERSUS SOIL DISPLACEMENT 

The difficulties inherent in the ELMS are illustrated in Figure 1. In the 
first place the choice of independent variable is not appropriate. To make the 
relation correct dimensionally, the slip, which is measured in inches, would 
have to be divided by some other length -- presumably the length of the test 
section originally used to obtain the empirical relation. Since no information 
is available on this point curve III on Figure 1 cannot be used for general 
analytical work. The second difficulty with the ELMS is that the friction 
coefficient T/a approaches zero for very small or very large displacements. 
When an attempt is made to use curve III in a graphical analysis of recipro- 
cating motion, it is found that for very short or very long strokes the drag 
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' tosELM? J ,‘ l * * * S " ot f cons ‘ dered reasonable. The third difficulty in using 
^ r f “ C n r r V reciprocating motions is that the friction coeffi- 
present- but fti n f "egative slippage. Some type of hysteresis must be 

o.Ven ’ b h f f he hyster esis loop cannot be deduced from the data 

motion th So!'l°L e „a f , 0 J , , t t inin ^ “ analy * ical <« the extension-contraction 

motion Soil Model II has been considered. This model is represented on 

Figm-e 1 by the straight line II through the origin, with the same slope as the 

ulate tbin° P th lmeUI ,- AUh ° Ugh tMs m ° del is ePPsiderably simpler to maSp- 
ulate than the original model, it also gives unrealistic results, and is not ? 
considered satisfactory. 

Since any reciprocating device disturbs the soil and can be considered to 
operate m something resembling loose sand after the first stroke, a third soil 
h^ef ' liri 6611 dev eloped based on the observed behavior of loose sand 
( . 0). These data presented m Figure 1 indicate that the friction coeffi- 

cient increases very rapidly with slippage up to a constant value which is 
independent of further displacement. This behavior suggests that the most 
ogical soil model would simply employ a constant friction coefficient, which 
s represented by line I in Figure 1. This soil model allows the desired per- 

b . r ““ C ® relations to be obtained in simple closed-form expressions, and has 
been used m all subsequent work. 

soil characteristic second only in importance to friction coefficient is 

surface tf 6 th S * * ° r •?* f °° tprint P ressure developed as penetration normal to the 
thl P?lfc 011 18 “" ea8ed - The range of soil hardnesses specified in 

shows Aat'the s^ "h \ A11 three 8011 m ° dels USe this Elation, which 

Shows that the soils being considered develop footprint pressures anoroxi- 

mateiy proportional to penetration depth. The proportionality constant, 

ver, varies by more than two orders of magnitude for different soils. 

l ij \ d ® fines f °°tprint pressure while penetration is increasing or 

held constant, but not while decreasing. This is a disadvantage, but not a 

bellows tvn * Te . Sei ? stud y* The distinguishing characteristic of 

laterTFil yP 7 QT 18 thelr low f °° tprint Pressure, which will be shown 

two-tenth.; J ] 15 * 6XpeCted to exceed 30 P° un ds per square foot or about 

f a P ° Und PCr square inch on the moon. Inspection of Figure 2 

shows that except in the case of the "extreme soil" the resulting penetration 

thereThe o! . n * gl ** ible and can be ignored completely. For the extreme soil, 
expected to t?n r ! 10 r may I appre u ciable ’ the an gle of internal friction is 
oUd tinrl l q •" th f 18 the soil behaves more like a liquid than a 

buoyant fore. SfTT varies directly with depth is equivalent to a 

as Jfunrtir, f the footprint pressure for deep penetrations will be considered 
etion of penetration depth only regardless of direction of motion. 
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PENETRATION DEPTH (INCHES) 

FIGURE 2. FOOTPRINT PRESSURE VERSUS PENETRATION 


EXTENSION- CONTRACTION TYPE OF VEHICLE 


The extension-contraction type of vehicle is a bellows vehicle that is pro- 
vided with a means of alternately increasing and decreasing its length, and 
derives its propulsion from this motion. The propulsion comes from a differ- 
ence in friction coefficient of the advancing and the retreating portion of the 
vehicle which is obtained by some sort of asymmetry in the construction of its 
lower surface. The primary reason for considering such a concept is the 
possibility of achieving a low footprint pressure device with a very simple 
mechanization, and with no exposed moving parts. The means of achieving 
the fore and aft differential in friction coefficient is discussed under Soil 
Considerations . 


Soil Considerations 

Producing the friction coefficient differential . - The most obvious way of 
attempting to cause the fore and aft differential in friction coefficient required 
for extension-contraction propulsion is to put cleats which slant in an aft 
direction on the bottom surface of the vehicle as in Figure 3(a). If the vehicle 
has a typical planform loading the result is likely to be as illustrated; both 
foward and aft moving cleats are fully imbedded in the soil. In a frictional 
material, this will result in the soil shearing along a plane at the lower edge 
of the cleats, and the resulting friction will be independent of the direction of 
motion or the inclination of the cleats. If, however, the cleats are built at a 
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shallow angle, and the vehicle is lightly loaded, the situation depicted in 
Figure 3(b) occurs. In this case, the forward- moving cleats ride on top of the 
soil and a friction coefficient differential is produced. This phenomenon will 
be referred to as planing, and the "cleats" will be referred to as planing sur- 
faces. The conditions for which planing will occur are discussed below, and 
it is shown that the primary requisite is a low vehicle planform loading. The 
disadvantage of this concept is its lack of versatility. It appears that a design 
can be conceived for almost any type of soil, but a given design will work in 
only a narrow range of soils. 

Another alternative is to hinge the cleats as in Figure 3(c). If the sur- 
face of the folded cleat is a normal rough surface, the friction coefficient of 
the forward- moving portion will be the tangent of the sand friction angle. 
Bekker^ indicates that the cleats on the aft moving portion may increase the 
friction by about 20 percent. For lightly loaded configurations, the weight of 
sand carried between the extended cleats is significant and the cleat augmen- 
tation may be greater than 20 percent. There is a possibility that the exposed 
surface of the folded cleats may be coated with a low-friction material such as 
tetrafluor ethylene. Tests of this material on sand made by the contractor 
showed a friction coefficient of about 0. 25. The folding cleats, shown in 
Figure 3, have the disadvantage of requiring exposed moving parts which is 
contrary to one of the primary advantages of the worm concept. 

Although a limited amount of analysis in this section will utilize the trac- 
tive effort as given by the ELMS model, most of the results are based upon a 
friction coefficient which is independent of the soil displacement as discussed 
in the section on Lunar Soil. The preceeding discussion applies in either case 
where, in considering the ELMS model, the term "friction coefficient" refers 
to the magnitude of tan <p. 

Theory of planing. - Using the Coulomb Earth Pressure theory as given 
by Jumikis,^ the geometry of a planing surface is shown in Figure 4. The 
surface is considered two-dimensional, with an angle of attack, ot, and is 
moving to the right. The internal friction angle of the sand is <p, and the 
friction angle of the sand against the surface is <pp The surface is supporting 
a running load of m p pounds mass per unit width (normal to the paper). It is 
desired to find the length, £, which will be in contact with the soil, and the 
horizontal force which will be required to move the planing surface. 

It is assumed that the triangular section of soil ABC moves with the 
planing surface. Shearing takes place along the plane AC, but the soil is just 
at the point of shearing on the plane AB. 


If the soil has a mass of m s pounds per unit width, the total weight being 
moved is (m p + m s )g/(g e ) 

per unit width, and the resistance is 
r - (mp + m s )g/( ge ) tan <p pounds per unit width. The effective friction 
coefficient therefore becomes: 

g r / m \ 

= — — = (1 + — - ) tan <P = tan (cp + a ) (1) 

6 gm p \ m p / 1 
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cos a + r 


sin a 


sin a 






J is the dimensionless planing parameter which relates the critical load on 
the planing surface to the soil density and the geometry of the planing surface. 
Since the critical load varies inversely with J, it is desirable to have the value 
of J as small as possible. To facilitate the use of Equation (4) in the design, 
the value of the dimensionless ratio J 2 y Q / rrtp is plotted as a function of 
the angle of attack O', and the friction angles <p and <p j in Figure 5. 

The curves of Figure 5 show that, if the planing surface is rough enough 
to force shearing to take place in the soil (<p j =< p), the value of J can never be 
less than about 4. However if the planing surface is very smooth or has a 
low-friction coating (<pi<( p ), the value of J can be considerably decreased, 
allowing higher loadings to be used. 


If the lower surface of a vehicle has a number of planing surfaces 
arranged in tandem, certain geometric constraints must be met for the 
planing analysis to be valid. The arrangement shown in Figure 6 is the 
maximum length of planing surface that can be effective on a vehicle 
surface. The problem is assumed to be two-dimensional. Increasing the 
length of the cleats will not increase the soil contact area because one cleat 
will shade the one behind to a degree determined by the soil friction angle, <p . 
The relationship between vehicle surface area and maximum cleat length is 


i 

v 


l cos a + 1 


sin Oi 
tan <p 


(5) 





(6) 


The ratio, R, approaches 1. 0 for small angle of attack, a, but is never less 
than 1. 0. 


The loading on the planing surfaces may now be referred to the planform 
loading of the vehicle. The vehicle planform loading, pp is defined as the 
total vehicle weight (in the local gravity) divided by the projected area of the 
vehicle bottom surface. Considering only the weight supported by one cleat, 
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M 

S 



(7) 


where small letters refer to quantities per unit of vehicle width. 
Equation (4) the allowable loading on the planing surface is 


m ~ 
P 



J 


lbm 


F rom 


16 


PLANING PARAMETER 



17 




FIGURE 6. PLANING SURFACES ON A VEHICLE BOTTOM 


per unit of vehicle width, 
surface is 


The corresponding pressure acting on the planing 


m 

= E • I_ 

& g. 


i y 


g_ 

g. 


( 8 ) 


This amount of vehicle loading will result in a vehicle planform -loading 

l y 


of 


P al P d ‘ 1 


RJ 


g_ 

g. 


(9) 


^r/ l0 " ble VeHCU P ' anf ° rm Planing 


be valid. 


To obtain a useful expression for allowable footprint pressure 
take a - 6 degrees, -Vs = 50 lbm/ft 3 = 0. 029 lbm/in. 3 , <p, = <p = 32 degrees 
ge/g = 6. From Figure 5, J = 4. 07. ° s 


R = cos 6° f 1 + tan 6 


fl + JSLiL\ . 1. 

\ tan 32° / 
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Then the maximum allowable vehicle footprint pre 


ssure is 


al 


nnsy 


(0.029 
(4. 07 


7*7 


i = 
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psi 


in. 


a V f hl *u le ls .. built wi . th a greater footprint pressure than this, 

a mound of sand ° r Wil1 assume . a positive angle of attack and 

T , .. . , will be built up under the vehicle as shown in Figure 7. 

e si ua ion depicted is a rigid object being towed through the sand. As the 
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FIGURE 7. LIGHTLY LOADED AND HEAVILY 
LOADED PLANING VEHICLES 

loading is increased above p a l, the vehicle assumes an angle of attack and a 
new planing situation develops with the entire bottom of the vehicle acting as 
a planing surface. The new allowable loading is based on the vehicle length, 
l, rather than the step length, l. 

£'y s g 

P al RJ ’ g (10) 

6 e 

Therefore, as the loading of a vehicle is increased above p a i, a pile of sand 
can be expected to begin to build up in front of the vehicle. 

Experimental data . - The experimental data shown in Figure 8 provides 
a rough verification of this simplified planing theory and indicates the con- 
sequence of exceeding the loading given by Equation (10). This test was con- 
ducted in dry sand with a density of approximately 100 lb/ft 3 , and a friction 
angle of 30 degrees. The model had a rough wood surface, and the geometry 
is shown in Figure 8. 


The allowable loading for planing was calculated by using J = 4. 0 as given 
y Figure 5. The differential in friction coefficient between the forward and aft 
directions disappears rapidly as p a l is exceeded, and is nearly nonexistent at 

a if£ aC v® ° f twice the allow able loading. Therefore, it may be concluded that 
although the expression for allowable loading, given by Equation (9), may not 
be entirely accurate, it is nevertheless a reasonable criterion for the loading 
ol an extension- contraction vehicle. ° 
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Extension-Contraction Vehicle Performance 


The analysis ol the mobility of the extension-contraction vehicles was 
made by equating the tractive forces for each end of the vehicle, i. e. , thrust 
equal to drag. \o dynamic forces were considered. The detailed analysis is 
given in Appendix C. The primary parameter used for comparing the pro- 
pu sive efficiency of vehicles is equivalent drag-to-weight ratio. This equiv- 
alent drag is not necessarily identifiable as an actual drag on the vehicle, but 
it is the drag which would consume the same power as the vehicle when 
applied at the average vehicle velocity. The equivalent drag-to-weight ratio 


■nn n c 


D 

e - powe r 

W W x velocity 


D 

e_ _ work 

W W x distance 


Three assumed soil characteristics were analyzed: 

C ase I . T/a constant - this is the classical friction coefficient where 
irictional force is proportional to normal force. This model is considered 

a good representation of loose sand if the vehicle stroke is large compared to 
uie cleat size* 

Cas^II. T/a proportional to soil displacement, j. This form is the first 
soil displ P a r ° Xlma t tl0n t0 the Engineerin g Lunar Mod el Surface (ELMS) for small 


soil displacement. 

Case III , r/a = C ( e “ 


0.1 


(elmst: 

1. 32 


J 


-0 


•* j ) 


■— \ . ^ u. i i d x iviuaei ^uriaci 

From Appendix B, the nominal value of C for a smooth surface is 


Engineering Lunar Model Surface 

m ~ r r _ 


stan/f HrtT 1 Pari ®° n . ° f ^ese three soil models was made by making the con- 
pi i <c , f coe t flcl ent for Case I equal to the peak value of T/a for the 
ELMS model; and the constant of proportionality for Case II, that which 

“f ' h u e lf Sl r ° f the ELM , S m ° del at the In each 

nominal multiplymg constant for the aft-direction was taken as 1. 2 times that 
e forward direction, following the reasoning given on page 14. 


The reason for studying Case II is that it gives information about the 
ELMS model since it is the first order approximation to the ELMS model for 
a small stroke. A complete solution to Case III was not obtained. 
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Power requirements for the three soil models. - The work required to 
extend a vehicle segment which has one end fixed to the ground is shown in 
Figure 9. This is obtained by integration of force times distance over the 
vehicle surface at any instant, and then integration over the stroke length as 
shown in Appendix C. Case II and Case III are parabolic with zero slope at 
the origin, and as expected, are equivalent for small stroke lengths. The 
friction coefficient assumed for Case I is equal to peak value of T/ct for the 
ELMS model. 



LENGTH OF STROKE (INCHES) F05403U 

FIGURE 9. WORK TO EXPAND A VEHICLE FIXED AT ONE END 


For Case I and Case II, the instantaneous stationary point, p, of a 
vehicle which is free to move remains fixed on the vehicle and on the ground 
during any half cycle of extension-contraction. The location of the point p is 
found by equating the forward force to the aft force on the vehicle. The 
motion of the vehicle midpoint (referred to here as the eg) is proportional to 
the distance from the point p to the eg. The work for the portion of the 
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vehicle on each side of the stationary point may be found as for Figure 9. 

•The equivalent drag-to-weight ratio is found by 

Work 
W W x s 

The assumption is made that each succeeding expansion or contraction is 
just like the first; that is, the soil displacement, j, is taken as starting from 
zero at the point at which the motion reverses direction. 

The result given in Appendix C is shown in Figure 10. The worst 
assumed condition for Case I is a friction coefficient equal to maximum 
T /a for the nominal ELMS model for the forward -moving portion, and Z0 per- 
cent greater than this for the aft -moving portion. The best assumed condition 
is for tetrafluor ethylene on the forward- moving portion. 

For Case III the stationary point p does not remain at a fixed point on the 
ground, but moves from its initial location. Then, since the soil shear 
stress is a function of the past history of the soil motion, the horizontal force 
is a complex function of the stroke and the location of the stationary point has 
not been solved. It appears that a rigorous solution requires a step-by-step 
numerical integration of the process. The initial location of the stationary 
point, however, is given as in Case II, and indications are that it does not 
depart greatly from this location up to stroke lengths of about 10 inches. 
Therefore, an approximate solution can be obtained by assuming that the sta- 
tionary point remains fixed throughout a stroke (half cycle). This is the 
result shown on Figure 10, and is accurate for small stroke lengths. It can 
be shown numerically that the initial motion of the stationary point is toward 
the rear of the vehicle (for an expansion stroke) and the vehicle will move 
further than predicted by the approximation. Therefore, at least the initial 
departure of the actual D e /W from the approximation shown in Figure 10 will 
be on the downside, i.e. , the actual efficiency will be better than predicted. 

Figure 10 shows that the drag-to-weight ratio using the ELMS model 
(Case III) may be made arbitrarily small by using either a very small or a 
very large stroke length. This occurs because the friction coefficient given 
by the ELMS model approaches zero for either small or large soil displace- 
ments. This causes serious doubts as to the validity or completeness of the 
soil model, at least when applied to the present situation. 

As discussed in the section on Lunar Soil, the results for Case I should 
be representative of the performance of the vehicle in loose sand, if the 
vehicle stroke is large. All further analysis and comparisons will be based 
upon this model. 

For conventional vehicles, the equivalent drag-to-weight ratio varies from 
a few percent (wheeled vehicle on smooth surface) to about 1. 0 for off-road 
operation in poor soil. For this study it will be assumed that the useful range 
of drag-to-weight ratios is below 0. 5. It should be noted that the installed 
power in any vehicle must be several times the propulsive power for operation 
on level ground to allow for acceleration, climb, transmission losses and 
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EQUIVALENT DRAG-TO-WEIGHT RATIO 






operation at less than full throttle. The general level of the drag-to-weight 
ratios shown on Figure 10 indicates that, at best, this type of device will be 
marginally acceptable, and that if the concept is to be successful, all avail- 
able means should be used to minimize power consumption. 

Slope capability . - A vehicle on a slope, 0, is shown in Figure 11. A 
constant friction coefficient, /i, as in Case I is assumed, and the subscripts 
f and a refer to the portion of the vehicle forward and aft of the stationary 
point. 



FIGURE 11. VEHICLE ON A SLOPE 
The equilibrium of forces parallel to the slope is 

M a w a = M f W f +Wtan9 (11) 


where, 


Substituting 


This reduces to 


or, 


where. 


W = W + W r 
a f 


L L 

W = -T-— W and W, = -J- W 
a L f L 


(/i a - tan 0) = j- (ji f + tan 0) 


L k = L. k. 
a a f f 


k = 


- tan 0 and k^ = + tan 0 


( 12 ) 


(13) 
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Equation (12) is identical to the expression given in Appendix C defining 
the relationship between L a and Lf. Other expressions from Appendix C can 
be shown to apply on a slope if the k a and kf, as given by Equation (13), are 
used. The resulting equivalent drag-to -weight ratio is 




(14) 


This is shown in Figure 12 for several assumed values of friction 
coefficient. The maximum slope a vehicle can climb is given by 


k = 
a 


k 


f 


M a - tan 8 = /Lt f + tan 8 

(tan 0) max = 1/2 > (15) 

This expression defines the asymptotes of the curves of Figure 12 since the 
equivalent drag -to -weight ratio approaches infinity at the maximum slope 
capability. This, of course, means that the forward travel goes to zero and 
the power required remains finite. 

From Appendix C, the force which must be exerted on the end of the 
vehicle to extend it against the frictional force is 


P 


1/2 W 


k k r 
a f 

k +k- 
a f 


(U ~ tan 0 ) (ix + tan 0) 

1/2 W _ * t 

+ M f 


( 16 ) 


Setting the derivative of this with respect to tan 8 equal to zero gives 

tan 8 = 1/2 (M a - M f ) (17) 

which happens to be the maximum slope the vehicle can climb. Putting this 
into Equation (16) gives the maximum value of the extension force: 

P max ff + M f > ( 18 ) 

Incremental e xtension. - In the preceding discussion, and in Appendix C, 
a vehicle with uniform strain throughout its length at all times has been 
assumed. The vehicle has been assumed to be constrained to this type of 
motion by some sort of internal mechanism, if necessary. A more realistic 
case might be a simple bellows with a spring force which is small compared 
with the friction force. In this case, each bellows convolution must have a 
restraining member to keep it from extending beyond its intended length. 

Also, in a vacuum, some device, such as a spring, must be provided to 
contract the vehicle. Now, if an extension stroke is begun with application of 
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opposed forces on each end of the vehicle, the bellows segments will pop out 
to their extended length, one at a time, as the magnitude of the force is 
increased. It may be shown that the distance traveled and the work required 
per stroke are the same for this case as for the uniform- strain case treated 
previously. The magnitude of the force required is not the same, however. 
For the uniform extension case, the extension force is a constant value of 


P 


1/2 W 


k k 


a 

TT+Tc 


a 


f 

f 


(19) 


throughout the stroke. For the incremental extension case, the extension (or 
contraction) force varies linearly from zero at the beginning of the stroke to 
twice the value given by Equation (19) at the end of the stroke. Because of the 
linear variation, the work required is the same for either case. 


CONSIDERATION FOR DESIGN APPLICATION 


The principal problem relating to the extension-contraction mode of 
operation is the design of a bellows lower surface which will achieve an 
adequate differential in shear between aft moving and the forward moving' 
portions. The general level of drag-to- weight ratios tend to be high and all 
means should be used to minimize power required. One such means is the 
use of a low friction type coating such as tetrafluor ethylene on the bellows 
lower surface during the forward stroke operation. A variable geometry 
profile such as a hinged or flexible cleat offers promise as it tends to retract 
during the forward stroke. A fixed geometry type bellows can be designed to 
plane and not bulldoze for any given soil but is not versatile for a range of 
soils. Footprint loading should be kept low. Configurations utilizing this 
mode of operation are restricted to forward motion only because of the 
ratchet-type operation inherent in the concept. 


RIB- WALKING TYPE OF VEHICLES 


The rib- walking concept consists of a circular motion of the vehicle ribs in 
which each rib or bellows convolution is 180 degrees out of phase with the 
adjacent one. This motion is analyzed below. In addition, two modifications 
to this motion are proposed in the following sections. The double-acting 
bellows is a rib-walking device with only three feet, and the traveling- wave 
has multiple moving ribs with less than 180 degrees phase angle between 
adjacent ribs. A mathematical model is given in Appendix D which may be 
used to describe any of these motions. 


Rib- walking 

If a rib- walking vehicle is operating on hard ground, and has no spring 
suspension, the motion of its center of gravity is as shown in Figure 13. 
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FIGURE 13. RIB -WALKING MOTION 


Each arc is the path described by one (e. g. , odd-numbered) set of ribs. 

At the tangent point of the arcs, that set of ribs leaves the ground, and the 
adjacent (e. g. , even-numbered) set takes over. 

It is assumed that the change in height shownin Figure 13 is accomplished 
without the expenditure of power, but that the kinetic energy of the system is 
lost each time the load is transferred to a new set of ribs. The tangential 
velocity is 


v = A 60 

but the average horizontal velocity of the vehicle is 



(19) 


( 20 ) 


The work required for each circular arc of Figure 13 is the kinetic 
energy lost at one impact 
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The equivalent drag-to-weight ratio is then 
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This result is shown in Figure 14. The efficiency, then, improves with 
increasing radius of motion, A. 



2 

AVERAGE HORIZONTAL VELOCITY (FT/ SEC) 
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FIGURE 14. EQUIVALENT DRAG-TO -WEIGHT RATIO 
OF A RIB -WALKING VEHICLE ON THE MOON 

A practical consideration of this motion is the roughness of the ride. 
The motion shown in Figure 13 involves a change in both vertical and hori- 
zontal velocity from zero to a maximum during each cycle. An additional 
consideration is that the vehicle will begin to jump at the velocity at which 
centrifugal force exceeds the vehicle weight and the analytical expressions 
will no longer apply. This velocity is given by 




and is plotted in Figure 15. Both Figures 14 and 15 indicate that this motion 
is limited to low speeds and that a large radius is desired. 

The effectiveness of the rib-walking motion may be improved if a non- 
circular path is used for the travel of the feet. However, this considerably 
complicates the mechanism and also reduces the clearance for the forward 
moving feet. 



RADIUS OF MOTION, A (FEET) F05408U 

FIGURE 15. RIB -WALKING VELOCITY LIMITATION 
IMPOSED BY JUMPING 


Double-acting Bellows 

The analysis of rib walking showed that performance improves with in- 
creasing radius of motion, and is independent of the number of ribs. This, 
then, led to the simplification shown in Figure 16, which is a mechanization 
which approximately simulates rib walking with only three ribs. This device 
is called a double-acting bellows, and consists simply of a crank at the center 
of the vehicle which moves a member that is attached to the two ends. The 
two ends move in phase with each other , and the center moves 180 degrees out of 
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phase with them. In the position shown in Figure 16, the center is moving 
forward, and the ends are approximately stationary. After a rotation of 
180 degrees, the vehicle will be curved concave upward, and the ends will be 
moving forward. 



FIGURE 16. DOUBLE-ACTING BELLOWS 

The exact form of this motion depends upon actual construction details, 
such as the bellows stiffness and the type of member used to attach the crank 
to the ends. 

In soft soil, the forward-moving portion of the vehicle will not be lifted 
completely clear of the ground, and will be slid along to a degree dependent 
on the degree of vehicle curvature and the hardness of the soil. In this 
respect it is related to an extension -contraction vehicle that partially lifts its 
forward-moving part. The mobility of this device has not been analyzed, but 
it is felt that its performance falls somewhere within the spectrum of the 
vehicles that have been analyzed. 


Considerations for Design Applications 

The rib-walking motion where adjacent ribs are 180 degrees out of phase 
has relatively low drag -to -weight ratios at very low forward velocities but not 
at modest velocities except where the path radius of the walking motion is large 
(which in turn increases the roughness of the ride due to high vertical dis- 
placements). The double-acting bellows concept offers a simplified approach 
for improving the performance of a bellows -type vehicle by achieving a large 
effective radius without a correspondingly high vertical displacement. 
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TRAVELING WAVE VEHICLE* 


For many years Aeronutronic has conducted studies of improved concepts 
for off-road vehicles based on analyses of mobility in nature. One of the 
first things learned in a study of animals is that most of them - in particular 
the ones which move close to the ground - tend to be long and slender. The 
extreme examples of the slenderizing trend are the snakes and various 
species of worms. These animals take advantage of their elongated form to 
spread their weight over a large area, and the resultant low ground pressure 
allows them to move over very soft ground. 

The methods of propulsion used by worm-like creatures vary consider- 
ably, but can be grouped generally into two classes. One class, represented 
by animals such as the earthworm, the measuring worm, and the leech, 
utilize an intermittent motion in which one part of the body and then another 
is moved forward. The other class, which includes the snakes, centipedes, 
millipedes, and most caterpillars use one portion of the body after another 
for propulsion, but the body as a whole travels forward with a smooth, 
continuous motion. 

Further study of the second class of animals just mentioned reveals 
that the propulsive movements in every case take the form of waves of 
motion which pass along the body, with several complete wavelengths over 
the body length. It is interesting to observe that in some animals these 
waves progress from head to tail and in some from tail to head. Only those 
animals in which the waves move from head to tail have the power of swim- 
ming, either through a fluid or between obstacles, with a snake-like motion. 
The effect of this traveling wave motion can be seen clearly when a snake is 
swimming on the surface of water as shown in Figure 17. 

One animal chosen for study by Aeronutronic was the centipede. Analysis 
of the centipede showed that its feet move in a manner which forms a wave 
pattern moving from head to tail, and that in doing this, they move approxi- 
mately in circles with nearly constant angular velocity and phase relationship. 
This observation prompted the construction of the mechanical model shown 
in Figure 18. This traveling wave model is made up of a series of cranks, 
rotating synchronously, with a 15 degree phase difference between cranks. 

This device can be used to simulate waves that either move aft as on a centi- 
pede, or forward as on a millipede, and has been very helpful in analyzing 
the possible motions. This previous work has been adapted to this study. 


Hard Surface Performance 

One way of reducing the losses of the rib -walking device on a hard 
surface is to make the phase angle between adjacent ribs less than 
180 degrees. The resulting motion for a phase angle of 45 degrees is 
shown in Figure 19, and is called a traveling wave. 

*See footnote on Page 1. 
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F03392 U 


FIGURE 17. SNAKE SWIMMING 


F03397 U 


FIGURE 18. TRAVELING WAVE MODEL 





The energy loss of this device on a hard surface may be seen from 
examination of its motion as shown in Figure 20. The phase angle between 
adjacent cranks is Ip. There is no discontinuity in the horizontal velocity at 
the point of transfer of the load from one crank to the next and only the kinetic 
energy due to the vertical component of velocity is assumed to be lost. 



FIGURE 20. MOTION OF A RIGID TRAVELING -WAVE DEVICE 

ON A HARD SURFACE 


The work required per arc of travel is kinetic energy lost at one transfer, 


Work = \ ~ (A CO sin J j j . 


(24) 


The average horizontal velocity is 
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and the distance traveled on each crank is 
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2 A sin 



Therefore, the equivalent drag-to -weight ratio is 
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For 0 = 180 degrees, this reduces to Equation (22), which was given for 

a rib-walking vehicle. It is shown in Figure 21 for moon gravity and a crank 
radius of one foot. 



VELOCITY (FT/SEC) F05413 U 


FIGURE 21. EQUIVALENT DR AG- TO- WEIGHT RATIO 
OF A TRAVELING WAVE VEHICLE ON 
A HARD SURFACE ON THE MOON 
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* The additional cranks in the wavelength are seen to greatly reduce the 
power losses as well as smooth the ride. Theoretically, this vehicle will 
begin to jump at the same velocity as the rib-walking device as shown in 
Figure 15. However, for a large number of cranks (0 small), the amplitude 
of the vertical motion becomes small, and small deflections of the vehicle and 
soil will serve to reduce the tendency to jump and to further smooth the ride. 


Soft Surface Performance 

It is seen that a wave as shown in Figure 19 moves toward the rear of 
the vehicle as the vehicle moves forward. It is possible to make the wave 
move in the same direction as the vehicle motion (at least on a hard surface) 
by turning the device in Figure 19 upside down. It is felt, for most applica- 
tions, that it is advantageous for the wave to move aft, since this minimizes 
the wave velocity with respect to the ground and thereby reduces the vehicle 
sensitivity to ground surface irregularities. Furthermore, in very soft 
soil, or even in a liquid, the vehicle is propelled by the swimming action of 
the rearward -moving wave. An analysis of this type of vehicle has been 
made for certain soil characteristics, and the results of this are given 
below. The details of this analysis are given in Appendix E. 

The analysis is based on the concept of Figure 19, but many cranks are 
assumed, giving a continuous, smooth surface, and the radius is assumed 
small compared to the wavelength, resulting in a sinusoidal wave shape. 
Therefore, kinetic energy losses of the type shown in Figure 21 need not be 
considered. 

The coefficient of friction between the soil and the wave surface is 
assumed to be a constant, fJL, which is independent of soil displacement. The 
normal force on the wave surface is assumed to be directly proportional to 
the depth below the original soil surface, regardless of the direction of the 
motion of the wave surface. The results of the analysis, based on these 
assumptions, are given in Figures 22 to 25. The results are normalized by 
the friction coefficient, jU, where possible, so that the results may be applied 
for any friction coefficient. All forces not acting directly on the wave surface 
are lumped as "external drag, " and may be actual drag forces or the com- 
ponent of the gravity force when on a slope, or a combination of these. 

Figure 22 shows the depth to which a wave will sink into the soil as a 
combined function of the planform loading, amplitude and soil pressure 
gradient, K<p. The performance in the following figures is given as a func- 
tion of this penetration parameter, z m /A. The vehicle velocity ratio is 
shown in Figure 23. This is the ratio of the vehicle velocity under the 
given conditions to the theoretical velocity on a hard surface. For reference, 
note that z m / A = 2 corresponds to a wave surface with its top just even with 
the soil surface. The equivalent drag -to-weight ratio is shown in Figure 24. 
The fact that D e /W is negative in some regions is due to the vehicle's ability 
to recover energy on a down-slope. Figure 25 shows the maximum slope 
capability of the vehicle. It is of interest to note that a vehicle that is sunk 
into the soil to the midpoint of the wave (where z m = A) can climb a slope 
equal to the friction coefficient. 
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Considerations for Design Applications 

The traveling wave motion can provide relatively low drag -to -weight ratios 
at all forward velocities. This ratio decreases significantly with decrease in 
phase angle between adjacent bellows segment particularly for high forward 
velocities. Decrease in phase angle also smooths the ride. It is advantageous 
to have the traveling wave move aft since this not only reduces the vehicle 
sensitivity to ground surface irregularities but also provides locomotion over 
and through very low-strength surface materials even when completely sub- 
merged although at a decrease in propulsive efficiency. 
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BELLOWS ANALYSIS 


Two major classes of bellows are of interest in this study. In one class, 
which will be referred to as corrugated bellows, the loads are taken by bend- 
ing stresses in the bellows skin. In the other class, which will be referred to 
as membrane bellows, the skin is thin and flexible, but reinforced by fibers or 
rings. In the latter class, most of the loads are taken by the membrane 
stresses which are produced by the internal pressure. In the following para- 
graphs some of the design equations for the two classes of bellows are 
summarized. 


CORRUGATED BELLOWS 


Corrugated bellows are made either from annular discs welded together 
along their inner and outer edges, or are formed by one of several methods 
from tubular blanks. The former process allows close spacing of the con- 
volutions and, therefore, a low spring rate. However, the resulting bellows 
is not well adapted to high pressure usage. The latter process produces a 
s moo thly -cu rved convolution shape which avoids stress concentrations. This 
shape is close to a shape which can take pressure stresses in tension rather 
than in bending. The bellows can adjust itself to high pressures by small 
changes in contour. It is, therefore, standard practice to ignore the pressure 
stresses in formed bellows of standard proportions. 

In the case of welded bellows, the shape is such that the pressure 
stresses are not significantly relieved by small deformations, and these 
stresses should be added to the stresses caused by deflection to obtain the 
total stress. 


Bellows Analysis 

Except for the difference noted and the changes in material properties 
caused by welding, the welded and formed bellows can be treated by the same 
methods. In this study the two-dimensional equations proposed by Den Hartog 
based on a suggestion by Timoschenko, will be used. These equations are 
developed and compared with test results in Reference 5. Although the 
equations apply exactly only to small extensions of flat-plate bellows whose 
depth is small compared to the bellows radius, they give good results on any 
bellows whose walls do not slope more than about 15 degrees and where the 
bellows depth is less than 15 percent of the bellows radius. This range covers 
all the applications of interest in this study. 

The following relations are for circular bellows of mean diameter D, con- 
volution depth d, and material thickness t. They apply to one convolution of a 
formed bellows or to a pair of discs in a welded bellows. 
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Stress due to a deflection 6 per convolution, with modulus of elasticity E: 


3 E t 6 
2d 3 


(28) 


Stress due to internal pressure is: 


£ . 

b 7T 


(29) 


The spring rate is: 


P 

s 


6 


E D t' 


2 d" 


(30) 


For a complete bellows of N convolutions, the spring rate will be: 


R = 
s 


P 

s 

W5 


E D t 3 
2 N d 3 


and the bending stiffness is: 


M 

0 


R D 
s 


3 3 

ff E D t 

16 N d 3 


(31) 


(32) 


A bellows of length L, whose ends are restrained against rotation and 
translation will buckle laterally at a critical internal pressure: 


2 ir r 
_ s 

^buckle L 


(33) 


In the case of a pressurized welded bellows, when it is desired to obtain 
the maximum deflection without allowing the total stress ( |f e | + l^b I) to 
exceed an allowable value, F, the optimum geometry is given by: 

( 34 ) 
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For a bellows so designed, the maximum deflection per convolution is: 


6 

max 


F 2 t 
3 Pi I 


(35) 


At this extension, the stress due to deflection, Equation (28), equals the 
stress due to pressure, Equation (29), and the total spring load in the bellows 


n D d 

= 6 Pi 


(36) 


That is, the spring load is equal to that caused by the pressure acting 
over 1/6 of the bellows depth. In most cases this is negligible. 


Design Conditions 

The conditions under which a bellows would actually be used in a vehicle 
are described in the following paragraphs. 

The bellows will be made up originally to some length L Q and then com- 
pressed to a shorter length L c for shipment. On arrival it will be expanded to 
some new length L 2 and will then operate about a point L 3 with a stroke not 
exceeding 6 n from Equation (35). Only during this last phase will the full 
internal pressure be present, and since for optimum design the stress caused 
by pressure equals that due to elongation, the latter will be limited to half the 
fatigue limit or about 12, 500 psi. This is based upon the assumptions that 
( 1 ) the pressure stress is considered, and is added to the deflection stress, 
and ( 2 ) the internal pressure increases and decreases simultaneously with 
bellows extension and contraction so that the same fatigue limit applies to 
both pressure stress and deflection stress. Both of these assumptions are 
approximately valid for an extension-contraction vehicle made of welded 
metal bellows. The operating points discussed are shown plotted approxi- 
mately to scale on a stress strain diagram for stainless steel in Figure 33 

Since the operating point L .3 in Figure 26 can be set at will, it will be 
located for the purpose of this study so that the maximum operating length 
(L3 + 6N) will fall within the linear operating range of the bellows, that is 

that the elongation per convolution be not greater than 50 percent of the ring 
depth. This range of travel is shown schematically in Figure 27. 

The closed length is equal to 2Nt, where N is the number of convolutions, 
and for optimum design, from Equation (34). 


t = 
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FIGURE 27. BELLOWS TRAVEL (SCHEMATIC) 
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Therefore, the closed length is given by 



and the maximum extended length by 



The effective thickness at maximum extended length is 


( 37 ) 


(38) 



The most probable shape for the bellows section is a flattened oval as 
shown in Figure 28, This provides maximum surface contact, and is reason- 
ably stable against overturning. For the purpose of this study it will be 
assumed that the section shape is such that the circumference is equal to 
3 times the flat width b. 
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FIGURE 28. BELLOWS SECTION 


For the shape shown in Figure 28, the most severe design condition is 
probably as indicated in Figure 29, where the section is supported at A and B 
and the horizontal span is subjected to bending from internal pressure, pay- 
load, and its own weight. 
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FIGURE 2 9. LOADING CONDITION 

To simplify the analysis it will be assumed that inflection points occur at 
A and B. Then, if conditions are uniform along the bellows and the useful 
load per unit length is called M u , the running load across the span AB per 
unit length is: 

r = p. + (-* + t . y ) i- lb/in. (40) 

p i \ b e / g e 

when p- is in psi, M is in pounds mass per inch along bellows axis, t g is in 
inches, 1 y is in pounds mass per cubic inch, and g/g e i- s r a-tio of local gravity 
to standard earth gravity (g/g e = 1/ 6 on moon). 

If the tension in member AB is neglected, the maximum bending moment 
at the center would be: 


M b 

max 



(41) 


and this is equal to the maximum resisting moment of a unit length of the 
bellows, which is: 


M b 


d t 

L = H e 

y d / 2" 


(42) 


Equating (41) and (42), and using the value of r from Equation (40) gives: 

2 t d 

= F — — (43) 

^ 6 

where t is given by Equation (39). 
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Solving for useful load, M 


M 


u 


T“ 


<*e ( 4 * e dF 

TV*-?- 


The empty weight is: 



M 

6 


3 b 

c 


( 44 ) 


(45) 


The useful load to gross weight ratio may be given as: 

M u 

Useful load = b 

Gross weight M ' M 

u + e 

b b 


(46) 


Substituting Equations (44) and (45) into this, using Equation (39) for t , 
and multiplying through by e 



this may be written as: 

Useful load = ^(b ) ^ ^ 

Gross weight f /d\ “ ~ b 


where: 


(47) 


f W ■ T 


F x 




(48) 


Design Charts 

In order to define the design limitations and the effect of parameter 
variations, design charts have been prepared. These charts apply specifi- 
cally to welded metal bellows, where pressure stresses are considered, and 
are based upon the equations presented previously. The stresses caused by 
pressure and deflection are assumed to vary together so that the same allow- 
able fatigue limit applies to both. 
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The derivation of the relations used for the charts and additional charts 
are given in Appendix F. Unless otherwise noted, the charts are based on 
stainless steel with a fatigue limit of 25, 000 psi, and a mass density of 
0.283 lbm/in.3 

Figure 30 shows the allowable elongation of vehicles made of stainless 
steel compared with fiberglass. Assumed material properties are shown. 

The allowable elongation is the property most effected by the difference in 
these materials. The parameters on the remainder of the charts are shown 
for stainless steel, but would not be much different if they were shown for 
fiberglass. 

Figure 31 defines a design condition which is useful in determining the 
limiting factors of lightweight design. The bellows depth -to -width ratio given 
here is for the vehicle which will structurally just contain the internal pres- 
sure; that is, it will withstand the bending stress in the bottom beam caused 
by internal pressure load, but not that caused by useful load weight or by its 
own weight. This condition is given by the solution of Equation (43) for 
y = M u = 0. It is clear that (under the design assumption that have been 

made) this vehicle will not quite work, even with no useful load. It can be 
shown, however, that only a small increase in d/b over the limiting case will 
allow the vehicle to structurally carry its own weight, and that a further fairly 
small increase will allow the vehicle to carry a very substantial payload. 

Figure 31 shows that a 60 percent useful load may be carried by a 400 inch 
wide vehicle with only slightly larger d/b than for the limiting case. The 
limiting condition is useful in being quite close to a useful design and the 
charts which follow are based upon this condition. 

Figures 32, 33, and 34 give the empty mass, the footprint pressure, 
and the material thickness of the limiting design condition. The value given 
for each of these parameters is the limiting value, and the values for a use- 
ful design are slightly higher. 

Figure 35 gives the ratio of empty mass to enclosed volume, which is a 
measure of the structural efficiency. The values for a cylinder and a sphere, 
shown here for reference , are much smaller than for the bellows devices. This 
is because the cylinder and sphere take the loads in tension, which is much 
more efficient than the bending mode in which the bellows operate. This 
comparison gives some idea of the penalty which is being paid for the mobility 
in the shelter designs. 

Figure 36 is based upon the planing criteria discussed in the mobility 
section. The components of the planing parameter on the abscissa are 
bellows material density, 7 ; soil density, 7 s ‘> R, the geometry factor from 
Equation (6), and J, given in Figure 5. The resultant pressure given is the 
maximum internal pressure to which a bellows vehicle can be designed if it is 
to meet the planing criterion. 
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FIGURE 31. ^ FOR THE VEHICLE WHICH WILL JUST CONTAIN 

THE INTERNAL PRESSURE 









400 



INTERNAL PRESSURE (PSI) 

F05422 U 


FIGURE 33. MINIMUM £MPTY VEHICLE FOOTPRINT PRESSURE 
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FIGURE 35. VEHICLE MASS TO VOLUME RATIO 
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MAXIMUM INTERNAL DESIGN PRESSURE 
IF PLANING IS TO OCCUR (PSI) 






Considerations for Design Applications 


In the design of corrugated bellows where loads are taken by bending 
stress in the skin, the pressure differential is a key factor in establishing 
geometry, skin gauge and structural weight, because of the local stresses 
introduced and requirements to prevent buckling. This is particularly 
important when the cross section of the bellows is noncircular. For a non- 
circular shape, the structural weight required to carry sizable payloads is 
very little more than that required for the empty structure. For convoluted 
bellows, the generally high stresses at the outer edge of the convolutions 
make difficult the attachment of any type foot of flexible cleat. 


MEMBRANE BELLOWS 


The membrane bellows refers to membrane -type bellows reinforced by 
rings or fibers. The membrane and the reinforcements are assumed to have 
negligible bending stiffnesses. 

A longitudinal section of such a bellows is shown in Figure 37. It is made 
up of a number of rings with radius r and spaced a distance "a" apart. Be- 
tween the rings is stretched a membrane which, because of the internal 
pressure, bulges out to approximately circular convolutions as shown. 



FIGURE 37. BELLOWS CONFIGURATION 
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• • If there is no external force acting on the bellows, it will extend to some 
length at which the internal forces are in equilibrium. To hold the bellows at 
any other length requires a force which depends on the internal pressure, 
geometry, and material properties. This spring force can be estimated as 
shown in the following paragraphs. 

Spring Force 

An enlarged section of the membrane is shown in Figure 38. The 
meridional section is assumed to be circular, with a radius Pj. The other 
principal radius at the attachment point to the ring is p 9 . 



The meridional tension force per unit width is q^ lb/ in. , and the corres 
ponding circumferential force From membrane theory we have: 


*1 + ^ 
Pi P 2 


(49) 
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The radial deflection of the membrane at the ring attachment is given by: 

6 = Et(<l 2 - ‘"’l) < 50 > 

where E is the elastic modulus, t is the thickness, and V is Poisson's ratio 
for the membrane. 

For most cases of practical interest, the ring will be designed so that its 
elongation is quite small compared to its radius. If this is the case, 

Equation (50) can be simplified by putting 6 = 0, giving: 

q 2 = ( 51 ) 


From geometry 


Pj = (r - r ) sec a 


and 


P 2 = r sec O! 

Substituting Equations (51), (52), and (53) into (49) gives: 


lr 


qi yqi 

-f 

r ) sec & r sec OL 
o 


P 


(52) 


(53) 


or 


p r(r - r ) sec a 

q l " r + V (r - r ) 

' o 


(54) 


The net longitudinal force at the ring section is the difference between the 
longitudinal pressure force and the longitudinal component of membrane ten- 
sion acting around the ring, 

^ 2 

P=p^7Tr - 2 tl r q^ cos a . 
or 


p. 77 r 

r i 



( 55 ) 
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• •The values given by this equation are plotted in Figure 39 for several values 
of V. 

In some of the bellows applications which are of interest, it is desirable 
to use an elastic membrane with inextensible reinforcements in only the 
meridional direction. For this case, 

q 2 = ° 


and 


q l ‘ Pi Pi 

= (r - r ) p. sec a 
o' *1 


(56) 


giving: 


P 


p. 7T r 


2 IT r(r- r ) p. 

o' *1 


or 


2 r 


p . IT r 


(57) 


This gives the same result as putting V = 0 in Equation (64). It is also 
the same as ignoring the curvature in the circumferential direction, and cal- 
culating the loads on the basis of two-dimensional, theory, as was done in the 
case of metallic bellows. 


The above relations can be compared with the method usually used to 
predict the spring force in an air spring which is: 


or 



(58) 


This relation is also shown plotted in Figure 39. All of the expressions give 
substantially the same result in the range which is of most interest for the 
present study, that is, for values of r Q /r greater than about 0.75. 
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Of the various methods for calculating spring force, the two-dimensional 
method is simple, and gives reasonably accurate results. Since it can be 
used directly to give local forces, or integrated to give moments, it will be 
used in the subsequent work. 

For use in design, it is desirable to convert the previous expressions 
which contain the unknown radius, r Q , into forms which are in terms of 
known quantities. For the cases which are of most interest, the meridional 
length ( {, in Figure 38) can be assumed to be a known constant. For this case, 
simple expressions can be obtained for the resulting loads. 

Using Equation (56), the longitudinal component of membrane tension 
loading per unit length around the ring is: 

q = qj cos a = Pj P i cos O' (59) 


From Figure 38, 


l - 2 OL P 1 


Putting Pj from Equation (60) into Equation *59) gives: 

n cos a 

1p~ ” 2 a 

The value of G! is obtained by noting in Figure 38 that 

a = 2 Pj sin a. 


Using Equation (60), this gives: 


l 


sin OL 

" oT 


(60) 


( 61 ) 


(62) 


(63) 


This relationship is shown in Figure 40, By using this relationship, 
Equation *611 may be graphically transformed into a function of a 111 , and the 
result is shown in Figure 41. 
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FIGURE 40. BELLOWS GEOMETRY RELATIONSHIP 


The spring rate can be obtained by differentiating Equations (61) and (63). 


i ( - OL sin a - cos O' \ 

~1P- ) 


cos a. - sm oc 


sin a. + 


cos a 


sin a 


- cos OL 


This is shown as a function of a/$ in Figure 42, the transformation being 
made by Equation (63). The force or spring rate per unit of circumference 
given in Figure 41 or 42 may be integrated around the ring to give total force 
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FIGURE 42. BELLOWS SPRING RATE 




,pr moment. If the rings are circular and parallel, the net axial load and 
spring force become: 


P = IT r p. - 2 IT r q 

= 2 1 r - if-] 


and 


( 65 ) 


= 2» r . p. = lit rp. . -La 
da da *i p. da 


( 66 ) 


where the functions q/U p^ and (1/pjHdq/da) can be read from the curves, and 
all other terms are known. 

The bending stiffness can be obtained from Equation (66) by integrating 
around the circumference as shown in Figure 43. 



FIGURE 43. GEOMETRY OF BELLOWS BENDING 


67 



For circular rings: 


da = 0 . r sin P 
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V 
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This can be compared with the section stiffness of a beam 
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(67) 


Buckling 


All bellows have a tendency to buckle laterally when pressurized inter- 
nally and restrained at the ends. The buckling load can be approximated by 
use of the Euler equation. 



( 68 ) 


In this equation, L is the column length between points of inflection, P can be 
replaced by its value from Equation (74) and El by its value from Equation (67) 
giving, 


r - cos a 
i OL 


TT 2 r 2 

2L 2 


sm a 
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sin a 


cos a 
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sin a 


- cos a 
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TT sin g; 
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COS 0! 
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r 

I 


1 

cos Ct 


a 


(69) 


It is interesting to note that according to this equation, the column length 
is independent of pressure since the stiffness increases with pressure. 


Shear Transfer to the Ground 

In an application of a membrane bellows, it is proposed to use the con- 
volutions as contact surfaces on the ground. For such an application, it is 
important to know how well the convolutions can transmit thrust to the sur- 
face, that is, what is the spring constant of a convolution in the horizontal 
direction. A simplified preliminary analysis of this characteristic is given 
in the following paragraphs. 

In Figure 44, points B and C represent the attachment points for the 
flexible membrane BDC, which carries an internal pressure p^ and forms a 
cylindrical surface with radius, p i. The surface is assumed to contact the 
ground at D. The spring force will be taken as the unbalance in membrane 
tensile loads at D when this point is moved horizontally through a small dis- 

nloromont r? <3 
1 
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FIGURE 44. GEOMETRY OF CONVOLUTION 
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The displacement of D to the right in Figure 44 causes the radius of the 
left-hand half of the cylinder to increase by d Pj and that of the right-hand 
half to decrease by d Pj. The unbalance of force per unit length normal to 
the paper is then the internal pressure times the net difference in radii, or: 

dP = p. • 2d Pj (70) 

From the geometry of Figure 44, ds equals l/(cos Ci/Z) times the change 
in length of BED. 


but: 


ds 


d(BED) 

a 

cos — 


(71) 


(BED) = 2 Pj sin ^ = 2 Pj sin 

( a - !> \ 
y since t ■ 4 y 

3X < BED > = 2 P] - -V cos TT7 + 2 5in 

1 ^ P I 


= 2 sin 


a 

T 


a 

ot cos 


(72) 


Using Equations (70), (71), and (72), the spring constant becomes: 


dP 

ds 


2 p^ cos y 

A- ( bed) 

d pj 


( 2 sin T 


- O' CO! 


(73) 


This result is shown in Figure 45 as a function of a I H , where again this trans 
formation is made by the use of Equation (63). This demonstrates the im- 
portance of keeping a/i large (close to 1.0) if the bellows is to be stiff in 
transmitting shear to the ground. 
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Consideration for Design Applications 


The membrane bellows where loads are taken by tensile stresses in the 
skin provide a minimum weight structure to withstand a given pressure dif- 
ferential. Buckling is prevented primarily by change in geometry and is little 
affected by pressure differential. Load can be transferred to the lunar sur- 
face through the covering membrane by means of the internal pressure. The 
membrane tends to conform to the surface being traversed and a very uniform 
load distribution should result. Attention is required in designing the skin- 
frame -linkage system to avoid excessive lost motion in transfer of shear to 
the ground. 



DESIGN STUDIES 


SCOPE 


Three preliminary studies have been made for the use of different bel- 
lows concepts as applied to several type missions in order to gain more 
insight regarding the desirability of using them in various practical situations. 
A different bellows concept has been studied for each of three type missions: 
lunar shelters which seldom have to be moved; a smail unmanned vehicle used 
to explore a limited area; and a manned, extended range roving vehicle. The 
requirements of these missions encompass a broad range of mobility condi- 
tions, and, while the concept selected for study in each mission may not nec- 
essarily be the most optimum, it is considered sufficiently applicable to 
indicate the potential value of that particular concept. 

The first study is for mobile shelters which would have the capability of 
moving from the landing site to the shelter site. Since ability to move is 
useful but not necessarily vital to the mission, the bellows concept used with 
shelters can have very low velocity, short range, low maneuverability, and a 
comparatively rough ride but weight of the mobility mechanism should be a 
small fraction of the overall shelter weight. The extension-contraction mode 
was selected as ucmg for thic ctudy. TVio at-nrUr is divided into 

two parts, one for a very large shelter and one for a small two-man shelter. 
Also briefly investigated is a novel concept of using the shelter as a very 
simple heat engine for propelling itself utilizing direct solar energy. 

The second study is for a small unmanned vehicle suitable for exploratory 
purposes. The bellows concept suitable for such use can have low velocity, 
limited range, a moderate amount of maneuverability and the ride need not 
be soft. The rib-walking mode utilizing the double-acting bellows principle 
was selected as appropriate for this study. While several different type 
double-acting bellows structures were considered and are described, the 
study was arbitrarily limited to a metal bellows configuration with the bel- 
lows structure not in contact with the lunar surface. 

The third study is for a manned, extended -range, roving vehicle. The 
bellows concept suitable for such use should have moderate velocity capa- 
bility, good propulsive efficiency, good maneuverability, ability to cope 
with a wide range of surface roughness and soil strength, and a soft ride. 

The traveling wave mode utilizing a membrane bellows structure was selected 
as appropriate to meet these comparatively restrictive requirements. 

All of the vehicles studied will require a steering and control system. 

The immediately following section describes such a system previously 
developed at Aeronutronic and considered applicable to control these vehicles 
in turning, climbing, and bridging. 
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STEERING AND CONTROL 


The wormlike creatures in nature control the direction of their motion by- 
bending their bodies. In many cases this control can be exerted in the verti- 
cal as well as the horizontal direction. An illustration of this is given in 
Figure 46, where a centipede is shown bridging across a gap in very rough 
simulated terrain. 



FIGURE 46. CENTIPEDE BRIDGING 


The centipede bridges by anticipating the changes ahead and programming 
its path to smooth out all but those of longest wavelength. Analysis and analog 
studies show that this technique has the capability of reducing average accel- 
erations on rough ground by an order of magnitude or more, for an average - 
sized vehicle. 

An attempt to mechanize some of these observed characteristics resulted 
in the Aeronutronic Bridging Off-Road Amphibian (BOA) concept. A powered 
model demonstrating this principle is shown in Figure 47. This vehicle uses 
a flexible, continuous track for propulsion, and is flexible and controllable 
both vertically and laterally. The vertical and lateral control is shown in 
Figure 48. The tension in the cables can be controlled remotely, thereby 
modifying the curvature of the vehicle. This is demonstrated in Figure 49 
with the model turning, and in Figure 50 with the model bridging. This model 
was able to negotiate trenches up to one -third, and walls of one -sixth of its 
length. In addition, the smoothing effect of the length was very noticeable 
even though the model had no flexibility in its suspension. 

The general method of steering and control which has been demonstrated 
with the Aeronutronic BOA model is applicable to the various configurations 
described in the following sections. 
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FIGURE 47. BOA MODEL 



CABLES FOR STEERING 


i 
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FIGURE 48. 


BOA STEERING MECHANISM 





FIGURE 50. BOA BRIDGING 
LUNAR SHELTER 


The lunar shelter would be used as living quarters by a group of astro- 
nauts for periods of up to one year. It would have limited mobility, that is, 
a range of a mile at a speed of 1/10 mile per hour. The sizes which are of 
interest range from a minimum two -man shelter up to the largest unit which 
could be carried as an Apollo payload (about 20 feet in diameter). 


Large Lunar Shelter 

Configuration requirements . - For this design, a large shelter which is 
normally stationary is desired. However, limited mobility - perhaps to 
change location once in its lifetime - is provided by the extension-contraction 
mode. An artist's sketch of the shelter is shown in Figure 51. 

A welded metal type bellows will be assumed and the design principles 
and nomenclature discussed in the Bellows Analysis section will be used. The 
following characteristics are desired 


Maximum width 

20 feet 

Maximum length 

60 feet 

Living pressure 

7 -i/2 psi 

Maximum speed 

0. 1 mph 
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1 mile 


Range 

Empty weight minimum 

A sketch of the shelter is shown in Figure 52. Three possible configura- 
tions for the corners are shown. A curved corner as in (a) avoids stress con- 
centrations, and the analysis used applies without modification if the radius 
of curvature is large. The welded corners in (b) and (c) are used for con- 
voluted bellows because of the difficulty of fabricating a curved corner. 
However, they involve unknown stress concentrations, and do not appear to 
have any advantage over (a) for welded bellows. 

Bellows design . - Initially the following characteristics are assumed: 
Semicircular cross section 


Material 


Stainless steel 


Fatigue limit 25, 000 psi 

Young's modulus 28.5 x 10^ psi 

Mass density 0.283 Ibm/in. ^ 

Gravity factor, g e /g t> 

For minimum empty weight, the minimum working pressure and mini- 
mum d/b (bellows depth/vehicle width) are desired. To a first approxima- 
tion, the working pressure can be one-fourth of the living pressure. This 
assumes that the "living" condition of the shelter is the unstrained position, 
and that while moving, it extends until pressure stresses are equal to bend- 
ing stresses. This is the condition resulting from optimum geometry as dis- 
cussed in the Bellows Analysis section. The allowable fatigue limit is 
assumed to be one-half of the allowable static stress. This is illustrated in 
Figure 53, where the stress is shown for one point on the surface of a bellows 
disc. The nominal length with working pressure is represented by A, and is 
seen to be one-half of the allowable fatigue stress. At this length, no stress 
due to extension exists, and the pressure may be raised a factor of four to 
point C, the allowable tensile stress of the material, if this pressure is not 
cyclic. Therefore, the design working stress must be greater than 
7.5/4 = 1. 87 psi so the desired living pressure can be contained. A working 
pressure of 2 psi will be used for the design. 

The design charts in Appendix F show that d/b (bellows depth to vehicle 
width ratio) of 0.04 will not contain a 2 psi internal pressure without exceedirg 
the fatigue limit. A d/b of 0.045 will contain this pressure and leave an 
allowance to carry its own weight and some payload. 
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(a) 


5.4 IN. 


(b) 

ALTERNATE CORNER CONFIGURATIONS 


I 


(c) 
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FIGURE 52. LARGE LUNAR SHELTER AND 

ALTERNATE CORNER CONFIGURATIONS 
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FIGURE 53. STRESS CONDITION OF BELLOWS 

The bellows configuration is now fully £££» 

r&^ttlnnnary of the assutne* 
and calculated values from Appendix F. 


Width, b 
Bellows depth, d 


Working pressure, p i 
Maximum extended length, L^ 
Material thickness, t 


240 inches 
10.8 inches 
0.045 
2 psi 

720 inches 
0. 0963 inch 
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Equivalent thickness, t 

0. 374 inch 

Number of convolutions, N 

128 

Empty mass, 

47, 000 lbm 

Useful load (mass), 

334, 000 lbm 

Loading density, M^/V 

37.7 lbm /ft 3 

Gross weight 

381, 000 lbm 

Nominal extended length 

675 inches 

Closed length (theoretical flat plate) 

24 inches 

Maximum stroke 

+45 inches 

Maximum footprint pressure 

0. 42 psi 

Maximum footprint when empty 

0. 052 psi 

Maximum pressure force less bellows (spring) 
force 

40, 000 pounds 

Bottom surface design. - If each convolution of the bellows ha s a fr.ru 
(planing surface) it might take the form shown in Figure 54 at the maximum 
extended length. 


FIGURE 54. 


H h— 


5.4 IN. 



BELLOWS AND FEET IN EXTENDED POSITION 
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unit length on these feet is 


At maximum gross weight, the load per 


m 


Gross mass _ g 


N 


381, 000 
= ijmrzzrm 

= 2.07 lb/ in. of vehicle width. 

... __4- +-V. ^ foot. it. 

If this load is assumed to be concentrated at me cem~ - - ” 

creates a moment of 


M, 


= 2.07^J^)= 5.6 in. -lb 


Which is transferred to a unit of running length of two bellows sheets. This 
creates a stress in the bellows of 


6M, 


f = 


(6)(5. 6) 


2? (2)(0. 0963) 


^ - 1810 psi 


comp a "^ to th. fatiw. limit of ZS, 000 p.i " ia ?. b ^ ha . 
absorbed in the conservator .of that • “ be^creaseli since they will 

flexible, the frmtion^coeff^cient^differeiU backward . 

ctre^ust be taken in the design, however, to ensure that the cavities 
between discs empty of sand on the forward stroke. 

Feet of this type may be initially stowed in the overlapped position as in 
Figure 55. 



FIGURE 55. FEET AND BELLOWS IN CLOSED POSITION 


Because of interference of the feet and their 

fold to the flat-plate thickness of the bellows material. As shown 
Figure 55, the vehicle closed length is 
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L 

c 



129 in. 


= 10.8 ft 


Mobility Analysis 


In the Mobility section, 
ing for planing is 


it was shown that the allowable planf orm load- 



Wo P . C Ca^c ^Uting 6 th e^e^r °f oot Wth^ “ all °)T able loadin g of 
present vehicles, gives gth ’ l, using the loading of the 

1 (desired f ° r Pining) = 420 inches (at maximum gross weight) 

= 52 inches (empty vehicle) 

shelter a factor of 100 with the 

that this applies onlv to the "tvr>‘ in e shelter empty. Note, however, 
that even tethe soil assumed The Z <* r “ular. y s = 50 lbm/f,3), ani 
this scale. The following analysi s D wi?l be hT?' not been verified on 
feet are properly designed to d a, >?e based on the assumption that the 
sary -efficient differential neces- 

As given in the Mobility section, the distance traveled in one cycle is 


2s = 2 L 


6 Vj* 

1 k a + k f 


2L 6. J- 

e lk 

a 

^7 


+ 1 


where, on a slope 




tan 0 


k f ~ + tan 0 

For S the ranle °" Fj gdre 56 for -veral assumed friction coefficient 

0.68 to 3 7ff feet on u”*? £rlctl °n coefficients, the travel per cycle is from 

level ground. This corresponds to 7800 to 1400 cycles 
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travel one mile. The time per cycle which will give the shelter a speed of 
0. 1 mph is 


T = (feet traveled per cycle)/0. 147 second 

This is shown on the right hand scale in Figure 56. 

The force required to extend the shelter against the soil frictional force 
on the maximum slope is 

/ M a + M f\ K 

^ -J( Gross mass) 

= (/Ll + /i )(0. 25)(0. 166>(381, 000) 

U I 

= 15,900 (/U + u) 

cl I 

This force is then between 15, 900 and 38, 100 pounds for the soil friction 
conditions assumed. It was shown previously that 40, 000 pounds of force 
was available at the working pressure of 2 psi. 


The equivalent drag -to -weight ratio was shown in the Mobility section to 
be 


D 

e 

w 


k k r 
a f 

k - k " 
a f 


Since the vehicle weight is now given, this may be converted either to 
power to travel at a given speed or work required to travel a given dis- 
tance. The power consumption is given by 

P = D v 
t e 

= (-f#)(M + M ) L JL v 

\W/ X e u' e g g 

= £(381,000) v(£) 

At 0. 1 mph (0. 147 fps), this becomes 

P t = 9330 ft - lb / sec 
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or 


k a k f 

HP = 17 £ — _ ' k horsepower 

a " K f 


This is shown in Figure 57. 


Small Lunar Shelter 

Configuration requirements. - The subject of this study is a minimum - 
size shelter for the protection of a two -man tean during an eight- to 
thirty-day stay on the moon; and includes the requirement that the empty 
shelter shall be able to move under its own power at low speeds over dis- 
tances up to about a mile. 

The design procedure used in this study is similar to that used in the 
previous section for a large lunar shelter. 

The nominal external dimensions for the shelter in a fully extended condi- 
tion are given in Figure 58, and a sketch is shown in Figure 59. It is de- 
signed for a normal internal living pressure of 5 psi, and is made up of a 
main shelter 7 feet in diameter by 12-1/2 feet long and an airlock of the same 
diameter, 5 feet long. 

As shown in Figure 58, the airlock and the shelter are of metallic 
bellows construction, and the entire assembly can be collapsed to a length of 
about 5 feet for transport. 

Bellows design. - For a final living pressure of 5 psi, the working 
pressure will be 1.25 psi. For this analysis it will be assumed that the 
bellows shape is equivalent to one of 60-inch width. Reading from Figure F-l 
of Appendix F, a d/b of about 0.037 will sustain a pressure of 1.25 psi so 
the bellows depth will be made. 

d = 0. 037 x 60 = 2. 23 inches 

Since the bellows will require welding, some of the material will be in 
the annealed condition, so an allowable stress of 25, 000 psi will be used. 

This stress level, together with the working pressure, defines the optimum 
geometry given in the Mobility section. 

i = ,4i 

t V Pj V 1-25 

The optimum thickness would therefore be 

t = = 0. 0158 inch 

141 
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FIGURE 58. SMALL LUNAR SHELTER LAYOUT 

In the subsequent work, the values of d and t will be rounded off to 

d = 2.25 inches ; 
t = 0. 016 inch 

If the maximum extension per plate is held at 25 percent of the depth, the 
bellows maximum length per plate will be 

0. 016 + 0. 25 x 2. 25 = 0. 579 

For the total active bellows length of 170 inches, the number of plates will be 


2 N 


170 
0. 579 


294 


The theoretical minimum closed length will be 

L - 294 x 0.016 = 4.7 inches 
c 

The ratio of maximum deflection to closed length is given by Equation (35) 
of the Bellows Analysis section 

F 2 (25, OOP) 2 

6 x 1. 25 x 28. 5 x 10 6 


6 

t 


3p.E 


2 . 92 



The total stroke will therefore be 

L = +2. 92 x 4. 7 = 20 inches 
s 

This is 44. 4 percent of the stroke calculated for the large lunar shelter, and 
the resulting travel per cycle is shown in Figure 60. This was obtained from 
Figure 56 by the required reduction of ordinate scale. 



FIGURE 60. TRAVEL PER CYCLE 


The mass of bellows material will be approximately 

M. = ft DL d y - x 84 x 4. 7 x 2. 25 x 0. 283 
l c 

= 790 lbm 

If the closing bulkheads are in the form of 2:1 ellipsoids of revolution, the 
maximum radius will be 

42 ^ 

r = = 84 inches 

m 21 
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If the maximum bulkhead stress is to be held to 25, 000 psi for an internal 
pressure of 5 psi, the bulkhead skin thickness must be 


t 


p.r 

l m 


F 


5 x 84 
2 x 25, 000 


0. 003 inch 


Since this is considerably less than would be suitable for handling loads, 
an effective thickness (based on projected area) of 0. 016 will be used, giving 
a mass for three bulkheads of 


M 2 = 3 x 0.785 x 84 2 x 0.016 x 0.283 = 7 5 lbm 


Doubling this mass to allow for reinforcing rings and attachment fittings, 
gives an order-of-magnitude estimate for the bare structural mass 

M 3 = 150 + 790 = 940 lbm 


To obtain traction in the extension- contraction mode, some type of foot 
structure will be needed. To a first approximation, the total foot area will 
be equal to the footprint area in the extended position, which is 

S = 60 x 170 = 10, 200 sq in. 

A lower limit for effective thickness for the feet would be about 0. 020, giving 
a mass for the feet of 


M . = 0.020 x 10,200 x 0.283 = 57. 5 lbm 
4 


Adding this to the structural weight of 940 lbm gives a figure of about 
1000 lbm for the total structure. 


At the extended length, the planform loading is 

1000 n nl / ik/* 2 

p - = 0. 016 lb/m. 

p (6) (10, 200) 

Then, according to the planing theory (see Equation (9)), for soil with density 
of 50 lbm/ft 3 , the minimum foot length for planing is 

H . ~ 1000 p = 16 inches 

min 1 p 

Building feet this large on this size vehicle appears to be a problem, but not 
necessarily an insurmountable one. In some types of soil, a foot on each 
bellows convolution would be adequate to produce a friction coefficient 
differential, and this could be accomplished quite easily. 
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Power analysis . - As in the case of travel per cycle, the power required 
to move the shelter can be obtained by reducing the values of Figure 57 in 
the ratio of masses, which is 1/381. The results are given in Figure 61 in 
terms of horsepower for a speed of 0. 1 mile per hour for various slopes and 
friction coefficients. 



SLOPE (PERCENT) 

FIGURE 61. POWER REQUIRED 
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As shown in Figure 61, the power required to move this shelter is quite 
low, but this power should be compared with that required to move it by other 
means. Since it is assumed that the shelter is empty and that some other 
method used to move its eventual contents, it is conceivable that the shelter 
could be dragged or rolled to its final destination. Assuming a friction 
coefficient of 0. 625, the power required to drag a mass of 1000 pounds on the 
moon at 0. 1 mile per hour would be 


HP 


1000 

6 


0. 625 


0 . 1 


375 


0. 028 
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If the shelter has a circular section and is rolled, the power requirement 
would be still less. 

The tentative conclusion from this brief study is that the extension- 
contraction mode is not the most economical method for moving a small 
shelter if other vehicles are available. 
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■ Solar -Mechanical Energy Transfer System for 

Extension -Contraction Vehicle 

A primary factor that makes the extension -contraction type of vehicle 
appear attractive is the possibility of utilizing the vehicle as a heat engine 
directly. This might be done if the vehicle can be made to selectively absorb 
or reject heat, resulting in a change in internal pressure due to a change in 
temperature of the gas inside. In a vacuum, the internal pressure can be 
used directly to extend the vehicle, but to contract the vehicle, some means 
such as a spring must be provided to supplement the low spring force of the 
beiiows structure. Also, a mechanism is probably required to hold the 
vehicle contracted while it warms up, and hold it extended while it cools off. 

The thermal property which is required is heat absorption when con- 
tracted and heat rejection when extended. Figure 62 shows the result of a 
preliminary thermal analysis of a configuration that has this property. The 
figure shows the equilibrium surface temperature as a function of the area 
ratio of two types of materials. 

The thermal properties of the surfaces are: 



Wflice SlilCd.Lt: pctxxxi. 


Hanovia gold | 6.0 

In the application of this concept it is expected that this area ratio is changed 
by the vehicle extension and contraction due to the change in geometry of the 
bellows. In this case, the two materials are Hanovia Gold 6518 paint (high 
temperature) and LMSC Silicate paint (low temperature). The assumptions 
made in the thermal analysis are: 

(1) Smooth cylinder 

(2) Longitudinal axis normal to the sun 

(3) The ends of the cylinder do not contribute to the heat transfer 

(4) An internal circulation system maintains a uniform temperature 
around the cylinder. 

Heat sources considered include the sun, lunar albedo (reflected solar), 
and lunar infrared. Diffuse reflection was assumed for the lunar noon case, 
and no reflection was assumed for the sunset case. Probably the most sig- 
nificant factor omitted from this analysis is the effect of the shape of the bel- 
lows corrugations. This could probably most easily be evaluated by exper- 
iment, but is beyond the scope of the present work. 

From Figure 62, under good conditions, maximum and minimum tem- 
peratures of 260°C and 90°C, respectively, are estimated. Regardless of 
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AREA OF GOLD + SILICATE/AREA OF GOLD 


FIGURE 62. SURFACE TEMPERATURE OF SMOOTH CYLINDER WITH 

GOLD AND SILICATE SURFACES 

the type of engine, its efficiency cannot exceed that of the Carnot cycle, 
which is 

„ _ T j " T o _ 260 - 90 "0. 313 

'c T. " 260 + 273 
1 

For a solar constant of 442 Btu/ft x hour the maximum possible power 
input per square foot of projected area would be 

P. = 0. 313 x 442 
in 

= 138 Btu/hr 


= 30 ft lb/sec 




If the sun is at or near the zenith the projected area is approximately 
equal to the footprint area, and the power required per unit area to move the 
vehicle at a speed v ft/ sec becomes: 


D 

P = -rjf? x p x v ft lb/ sec 
req W 


where p is the footprint pressure in pounds/ square foot. 

In an actual case, Pi n will be reduced by a number of factors, such as 
the actual cycle efficiency, the ratio of projected area to footprint area, 
mechanical losses, etc. If all these are lumped together in a general effi- 
ciency factor n, the resulting velocity becomes 


_ T| x 30 

(t) x p 


ft/ sec 


Values of v, for different values of T) and p and for an assumed value of 
(D e / W) equal to 0. 5 are plotted in Figure 63. The figure emphasizes the 
importance of low footprint pressure, and gives an idea of what overall 
efficiency must be achieved to give useful performance. 
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FIGURE 63. ESTIMATED PERFORMANCE OF EXTENSION- CONTRACTION 
BELLOWS VEHICLE UTILIZING SOLAR-MECHANICAL 
ENERGY TRANSFER SYSTEM 
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Design Evaluation 


An all-metal bellows design for a lunar shelter does not appear to present 
any unusual problems regardless of shelter size. The shelter can be designed 
with footprint loading sufficiently low to allow planing for small size shelter 
but not for the large shelter by an estimated one order of magnitude. It is not 
clear however whether this is really the problem it seems, since the planing 
theory has not been confirmed, and the true soil properties are not known. The 
extension-contraction mode is not a particularly economical method for pro- 
pulsion with regard to power consumption. The solar -mechanical energy 
transfer system however holds promise as a very inexpensive, lightweight 
means for supplying the mechanical energy for such operation. 


SMALL UNMANNED VEHICLE 

This section contains a preliminary design study of a small unmanned 
exploration vehicle utilizing the double-acting bellows concept, employing a 
circular convoluted stainless steel bellows. 


Configuration Requirements 

The following design requirements for the vehicle have been selected: 


Useful load, including power plant 

100 

lbm 

Average speed 

1 

mph 

Total range 

10 

miles 

Total life 

15 

days 

Diameter 

3 

feet 

Bellows working length 

9 

feet 


Possible Configurations 

The double-acting bellows principle selected for the small manned vehicle 
study can be utilized in a variety of configurations. Type of bellows structure 
and method of achieving roll stability are two variables. Several configura- 
tions are briefly described below. 

Flat-bottom welded metal bellows . - Sketch (a) of Figure 64 shows a 
double-acting bellows of welded construction. Since the welded bellows, 
unlike the convoluted bellows, can be made with a noncircular section, it is 
possible to use a flat bottom and eliminate the feet necessary with a circular 
bellows and thus obtain roll stability. The bottom surface could be coated 
with an elastic, abrasion-resistant coating to prevent damage and wear. 
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Practically the whole planform of the vehicle would be available as footprint 
area on soft surfaces. On hard surfaces, the contact area would be localized, 
thus minimizing slippage. 

Side-by-side circular bellows . - Sketch (b) in Figure 64 shows another 
arrangement for achieving roll stability. Here, two circular bellows are 
used side-by-side. This is probably the simplest way to create a small 
vehicle, since standard components could be used. The lower surface of the 
bellows would require protection, and, because of the possibility of high 
local loads, the concept probably would be applicable only to small vehicles. 
Subject to this limitation, the concept appears to be quite practicable. 

Two-lobed diaphragm bellows. - Sketch (c) in Figure 64 shows a possible 
application o7 the membrane bellows concept to the double-acting bellows con- 
figuration. Roll stability is achieved by the use of a two-lobed section, and 
since the entire structure is loaded entirely in tension, it should be very light. 
In a vehicle of this type there might be a problem in attachment of concen- 
trated loads, so this configuration might be suitable for transporting bulk 
loads such as liquid or gas. Since the structural weight of this configuration 
approaches very closely the minimum possible weight for a pressure vessel, 
it has possibilities as a mobile storage container in very large sizes. 

Circular bellows with feet, - Figure 65 shows an arrangement of cir- 
cular metal bellows with roll stability obtained by the use of feet attached to 
the end bulkheads and to a non-bellows midsection. The arrangement, while 
probably not the lightest structure, offers no unusual problems beyond the 
state-of-the-art of bellows technology. The circular shape bellows are not 
subjected to abrasion or unusual eccentric loads caused by direct contact with 
the lunar surface such as would be experienced by the three previously con- 
sidered arrangements. It was decided to limit to this arrangement the small 
unmanned vehicle study detailed in the following sections. 


Bellows Design 

The double-acting bellows utilizes a variation of the rib- walking motion, 
in which the path of travel of the feet is not necessarily circular. However, a 
preliminary choice of the stroke can be obtained from Figure 14 of the 
Mobility Studies section. This figure shows that in order to hold the effective 
drag- to- weight ratio in the ’’useful range” below 0. 5 for a speed of 1, 5 feet 
per second, the path radius should be at least 0. 5 foot. Based on this, the 
crank radius will be set at 7 inches. With this radius, the maximum travel 
per crank rotation, assuming no slippage, would be 14 inches. Actually there 
will always be some slip, and in soft soil this can become quite large. 
Therefore, an average travel of 7 inches per cycle will be used for this study. 
This means that for the specified total range of 10 miles, close to 
100, 000 cycles will be required. 

inn rwf cha . rt ln Fl g ur e 66 is based on test data and shows that for a life of 
, n ’ cycles, the completely reversed bending stress should not exceed 
, 0 psi. Since the pressure in the bellows will be constant during 
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operation, the local bending stresses due to pressure will be ignored. 

Because the hoop stress due to pressure is negligible in this vehicle, it 
also will be ignored. 

The configuration chosen for the design is shown in Figure 67. It is made 
of a short cylindrical center section which houses the payload, equipment, and 
power source, and two domed ends which are connected to the central section 
by lengths of circular bellows. Each of the three main components is sup- 
ported by feet which rest on the ground. The center section supports two 
transverse shafts which have cranks of 7-inch radius on each end. The shafts 
are driven in such a manner that all four cranks are parallel at all times. 

The end bulkheads are connected to the cranks by four parallel cables which 
normally are equal in length so that while displacements normal to the vehicle 
centerline are possible, the three sections must remain parallel to each other 
at all times. 

The motion of the vehicle is similar to that described under Double- 
Acting Bellows in the Mobility section, except that since four cables are 
used instead of two, the rotation of the end bulkheads is eliminated. This 
increases the efficiency, since sliding of the feet is reduced. The use of 
four cables also prevents instability since the additional cables force the 
bellows to buckle in a higher mode. This is illustrated in Figure 68, which 
shows the effect of various restraints on the buckling modes. 

In Figure 68, the lowest buckling mode of a fixed-ended bellows is shown 
at (a). For this mode, the equivalent pin-ended column length, L (the dis- 
tance between points of inflection), is one-half the distance between supports. 
If one of the ends is allowed to move laterally without rotation, as at (b), the 
equivalent column length is doubled and the bellows will buckle at one-fourth 
of the buckling load of (a). If, however, the median plane of the bellows is 
not allowed to rotate, as at (c), the equivalent column length again becomes 
one-half the overall distance between the end supports, and the buckling load 
is the same as at (a). This shows that when the three sections are held 
parallel, as at (c). Equation (33) for the critical buckling pressure is 
applicable. 

The next consideration concerning the bellows is the amount of strain to 
be absorbed. The two critical conditions are shown in Figure 69 (a) and (b). 

In (a), the crank is horizontal and the average strain in the bellows is 2 A/ L. 
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At (b) there is no longitudinal strain at the centerline, but the vertical deflec- 
tion, z, at the center causes the bellows to bend to a radius which reaches a 
minimum value of L 2 /24 z at each end. This curvature produces maximum 
strains in the bellows of 12 Dz/L 2 . If the strain caused by bending is to be 
held to the same value as the longitudinal strain at the centerline, the ver- 
tical deflection must not exceed the value given by 

12 Dz 2 A 


or 

A L 
z = 

6 D 

7 x 108 0 - . , 

- = 3. 5 inches 

6 x 36 

That is, if the bellows fatigue life is not to be shortened, the vertical 
deflection should not exceed half the horizontal deflection. In addition, some 
parts of the bellows (the top and bottom portions at each end) may receive 
more than one stress reversal per crank cycle, with a possible lowering of 
fatigue life. Because of the many other unknowns in the study, it does not 
appear worthwhile to go into further detail on this point at present. 


The minimum propulsive energy would be required if, on the forward 
stroke, the feet are slid along the ground, but with zero contact pressure; 
therefore the internal pressure and the crank location will be proportioned 
to accomplish this as nearly as possible. 

Based on previous work, it will be assumed that the total bellows mass 
is 100 pounds, giving a gross mass of 200 pounds, of which 133 pounds is 


supported at the center and 33 pounds at 
placed below the bellows centerline by a 
extremes of vertical cable displacement 
Figure 70. 



FIGURE 70. CABLE 


each end. If the crank axis is dis- 
distance y, the condition at the 
are as shown at (a) and (b) in 



DISPLACEMENT 
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Assuming that the mean bellows diameter is 35 inches, the cable load 
will be 


P = 1 X 35 2 p. = 960 p. lb 
4 *i 


and the two conditions for Figure 70 are: 

for (a): 960 p. x i^l = 21 


and for (b): 


2 x ,60 p. x ll±Xl = ill 

1 54 6 


Solving these equations gives: 


y = 2. 33 inches 
= 0. 07 psi 

A pressure of 0. 1 psi will be usea in tne renidiuiug w«^xx^. 

The choice of the actual bellows convolution shape is based on two 
requirements: the required deflection, and the necessity of preventing column 
or "squirm" instability. The equivalent fixed-ended column has a length of 
108 inches, and the required minimum total spring rate for the entire bellows 
is given by Equation (33). 


^buckle 


2^ (R s ) 

L 


or 


R 

s 


108 x 0. 1 
2 


1. 72 lb/in. 


This rate is so low that it would be exceeded by any useful bellows, so the 
only condition remaining is the required deflection. A bellows depth of 
1. 25 inches is within the range of good practice for a bellows of this diameter. 
The other bellows parameters are: the inside radius of the convolution to be 

5 times the skin thickness, and the slope of the bellows walls to be 15 degrees. 
A sketch of the resulting shape for a skin thickness of 0. 008 inch is given 
in Figure 7 1 . 
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From Figure 72, based on computer calculations, the maximum allowable 
deflection per convolution for the bellows of Figure 71 is 

6 = 0. 13 inch 

Since the bellows length per convolution is 0.825 inch, the number of 
convolutions for each bellows section, assuming a convoluted length of 
50 inches, is: 


N 


50 

0. 825 


^ 60 


and the allowable strain is 

N5 - 60 x 0. 13 = 7. 8 inches 



F05454 U 

FIGURE 71. BELLOWS CONVOLUTION SHAPE 


From Figure 73, the spring rate per convolution is 400 lb/in. giving 
a resultant linear strain rate for the whole bellows of 


R 


400 


3. 33 lb/in. 


BELLOWS DEPTH, d (INCHES) 


0 .05! 


0.1 
0.01 



- 


NOTE: DEFLECTION VALUES OBTAINED FROM THIS 

GRAPH ARE THOSE WHICH WILL PRODUCE A . 


MAXIMUM CONVOLUTION BENDING STRESS 
OF 30 000 PSI IN EITHER EXTENSION OR 
COMPRESSION OF THE BELLOWS 

- 

- 

- 

- 

1 1 1 1 1 1 1 1 1.JL 

i 111 i l l i i i 


0 . 10 

DEFLECTION PER CONVOLUTION, b (INCHES) 


1.0 


F05455 U 


FIGURE 72. MAXIMUM DEFLECTION PER CONVOLUTION 
VERSUS BELLOWS DEPTH 
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This indicates that the 0. 008 inch skin thickness could be reduced to 


t 



0. 0058 


without instability occurring. 

The bellows mass per convolution for a wall thickness of 0. 008 inch, from 
Figure 74, is 0.067 pound, giving a total mass for both bellows sections of: 


M, = 0. 067 x 120 = 80. 5 lbm 

bellows 


This appears to be consistent with the original estimate of 100 lbm for 
structural weight, especially w hen it is considered that tuc s Kin tuiCKiiess 
could be reduced considerably without incurring instability. 


The collapsed length of one convolution is about 0. 12 inch, giving a total 
collapsed length per bellows of 


L = 60 x 0. 12 = 7.2 inches 
c 


Since this is probably much less than half the length of the payload sec- 
tion, the latter will govern the collapsed length of the vehicle. 

The vehicle, shown in Figure 67, has a loaded mass of 200 lbm, giving 
it a footprint pressure on the moon of about 4. 2 pounds per square foot when 
all feet are on the ground. The power requirement will depend on the surface 
over which it travels, but at the target value of D e /W = 0. 5, the power con- 
sumed by the propulsive system at 1 mile/ hour would be 


P 


200 0.5 x 1.47 

~E~ X 550 


= 0. 0445 HP 
= 33 watts 


Steering of this vehicle would be done by differential adjustment of the 
longitudinal cables. The vehicle has the capability of moving backward as 
well as forward, and can be collapsed at any time to about one-third of its 
normal length. 


Design Evaluation 

The double-acting bellows configuration which was studied is considered 
to be quite suitable for meeting the specified mission objectives. Structural 
weight and power requirements are considered reasonable. The only problem 
apparent at this time is not in the bellows structure but in the design of the 
feet, which undoubtedly will require more study to derive a configuration 
suitable for use on the surfaces expected to be encountered. Because of its 
simplicity, absence of external mechanism, and ease of deployment, the con- 
figuration should be quite suitable for many unmanned missions. 
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MANNED EXTENDED-RANGE ROVING VEHICLE 


Configuration Requirements 


This section contains a preliminary design study of an extended range 
roving vehicle suitable for use by two men to carry out traverses of the lunar 
surface or perhaps to transport supplies. It utilizes the traveling wave mode 
with a membrane bellows structure. An artist's concept of this vehicle is 
shown in Figure 75. The vehicle is flexible and steerable both horizontally 
and vertically and can be collapsed for shipment or storage. It produces a 
smooth, constant-velocity motion, and appears capable of high speed. The 
configuration requirements selected for it are as follows: 


Useful load 
Average speed 
Minimum range 
Minimum life 
Diameter of crew cabin 
Length of crew cabin 


1000 lbm 

5 mph 
100 miles 

6 months 

6 to 10 feet 
10 feet 


Bellows ijesign 

The details of this vehicle are shown in Figure 76. It is made up of two 
coaxial bellows of which the inner one is pressurized to a living atmosphere, 
such as 5 psi, and forms the main structure of the vehicle. This inner 
bellows is reinforced by frames which determine its shape and support the 
payload and equipment. Also, each frame carries a transverse shaft with a 
crank on each end. These cranks are driven in synchronism and fixed phase 
relationship by a power transmission system, and each crank carries with it 
an outer frame to which is attached the outer bellows. 

The outer bellows is of flexible material and is pressurized to a compar- 
atively low pressure. The movement of the outer frames causes this bellows 
to deform in a manner which supports and propels the vehicle. Since the 
main structural member is a bellows, it has some flexibility and can be 
warped horizontally or vertically by cables running along the sides as in the 
Aeronutronic BOA model. 
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FIGURE 75. LUNAR SUPPLY VEHICLE 






Referring again to the vehicle configuration sketch in Figure 76 it can be 
seen that the shape of the lower surface is not the sine wave which might have 
been expected. Rather, it is a curtate cycloid with the lower portion more 
sharply curved than the upper portion. In addition, the lower portion 
skin is compressed while the upper portion is expanded. This is ^ ot the 
result of the particular linkage which is used, but is a fundamental char 
istic of this type of wave motion. 


In the configuration which is being studied at present, the skin of the 
outer bellows is being used to transfer the loads to the ground by means of 
the internal pressure of the bellows and the resulting membrane stresses. 
Since the skin is compressed longitudinally at the feet, it can billow outward 
at these points forming soft cushions to support and propel the vehicle. 
However, one problem which arises here is that these pillows are quite flex- 
ible in shear, and this causes some lost motion. The shear stiffness of the 
pillows is investigated in a preliminary way in the section on membrane 
bellows which indicates the nature of this problem. 


From the studies made on the traveling wave concept, it appears that the 
larger the number of cranks or the equivalent that are used per wavelengt , 
the smoother and more efficient the resulting motion will be. In the propose 
design, eight cranks are used per wavelength. At the present time it cannot 
be said with certainty that this is the optimum number. 


For an assumed effective drag-to-weight ratio of 0. 25, which should be 
quite feasible with such a vehicle, the propulsive power required on the moon 
would be: 


P 


3000 

E~ 


x 


0. 25 x 


5 

375 


1. 65 HP 
1.23 kw 


This power seems quite low and could be supplied by any one of a number 
of power sources. 

The footprint pressure of a vehicle of this type adjusts itself to the sur- 
face over which it is traveling. On hard ground the footprint pressure equals 
the pressure in the outer bellows, which in the configuration shown is about 
0. 35 pound per square inch. On soft soil the entire lower surface of the 
vehicle can rest on the surface, and the average pressure would drop to about 
0. 02 pound per square inch or about 3 pounds per square foot. As shown m 
the section on Mobility, the vehicle would still be able to make progress 
under these conditions, but the speed would be reduced by about 50 percent 
for the same power input. 
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Design Evaluation 


A traveling wave bellows design concept utilizing membrane construction 
appears quite attractive for a manned extended- range roving vehicle. Both 
structural weight and power requirements are considered quite low. The 
vehicle tends to mold itself to the terrain with a very low footprint pressure 
which should permit travel: (1) over a spiny type surface (such as a bed of 

nails), (2) over and through very soft soils, and (3) over rough terrain with- 
out vehicle lift-off and bounce under low gravity conditions. Its feasibility is 
dependent upon the availability of a suitable membrane covering material. 

This cover must withstand the lunar environment while still remaining flexible 
and abrasion resistant. The concept produces a smooth, constant- velocity 
motion and appears capable of moderate to high speed and efficiency as well 
as great versatility. 



OVERALL EVALUATION 


With any technical concept, the ultimate measure of practicability of the 
concept lies in its ability to perform some mission better than any other 
concept. The result depends upon the mission definition and its specific 
requirements as well as the measures of performance that are used to eval- 
uate and compare the systems. A detailed evaluation is beyond the scope of 
this study since the input data are not sufficiently well defined. However, a 
few observations will be made here and an attempt will be made to put the 
lunar worm vehicles into perspective with respect to more conventional 
vehicles. 


COMPARISON WITH OTHER VEHICLES 


Vehicle Types 

To realistically evaluate the new concepts which have been presented in 
this study, parametric comparisons will be made between those concepts and 
more conventional vehicle types. These types are listed and illustrated by 
diagrammatic sketches in Figure 77. 


Vehicle Parameters 

Vehicle reliability is of prime importance in a lunar vehicle. The degree 
of reliability of a vehicle mechanism is not as much of a characteristic of 
vehicle type as it is of the amount of detail design and development effort 
which have been expended on it. Reliability in the sense of probability of 
completing a mobility mission, however, is related to some of the basic 
capabilities of the vehicle. 

Since there is a good possibility that allowable soil bearing pressures 
may be low on the moon, it seems logical to classify vehicle types with regard 
to the degree to which their weight is distributed over the total planform 
envelope of the vehicle. In general, on soft or fluid surfaces it is necessary 
to spread the weight over as large an area as possible, and this is facilitated 
by the fact that such surfaces are usually reasonably smooth and level. With 
increasing roughness of the surface, the portion of the vehicle in contact 
with the surface necessarily decreases, but it still is an advantage to have a 
large contact area because of the averaging effects of widely spaced contact 
points. If the surface is rough but comparatively weak structurally, a dis- 
tinct possibility on the moon, many of the higher projections will be flattened 
down by the vehicle, particularly if the contact surface on the vehicle is 
unyielding. Although this process tends to smooth the ride it is undesirable 
since it consumes considerable power; therefore, in this it is also beneficial 
to increase the contact surface area. Only in the case of a hard smooth sur- 
face is it desirable to use small contact areas and the accompanying high 
footprint pressures. 
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WHEELED 





Based on the above reasoning, the ratio of footprint area to planform 
envelope is considered a basic parameter for comparison of vehicles. For 
any unprepared surface, it is desirable to have the footprint area as large as 
possible. 


In order to be able to travel over any type of surface on the moon, a 
parameter of major importance in a vehicle is the size of obstacle it can 
clear. Since obstacle clearance usually requires climbing, the vehicle auto- 
matically will be able to climb a slope equal to the obstacle height/ vehicle 
length, therefore, this latter ratio will be used to characterize the vehicle 
maneuverability. The various vehicle systems can be classified generally 
with respect to footprint area and obstacle clearance as shown in Figure 78, 
where each type of vehicle is shown occupying a shaded area on the chart. 

The values for the conventional vehicles are typical values from references 
such as the OTAC Tracked Versus Wheeled Vehicle report. ° The values 
for the worm type vehicles are estimated on the basis of the analysis and 
design studies in this report. The values obtained from Figure 78, together 
with the fineness ratios (length-to -width ratio) of the vehicles are listed in 
e o owing table and from design studies of propose'd lunar wheeled vehicles, 
References 8 and 9. 



FIGURE 78. VEHICLE OBSTACLE CLEARANCE CAPABILITY 
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Vehicle Type 


Wheeled 

0. 03 

- 0. 06 

0. 005 

- 0. 0i3 

3 

Tracked 

0. 1 

- 0. 15 

0. 15 

- 0. 25 

2 

Tandem Tracked 

0.07 

- 0. 12 

0. 6 

- 0. 8 

4 

Extension- Contraction 

0. 01 

- 0. 04 

0.8 

- 0. 9 

6 

Traveling Wave 

0. 13 

- 0. 17 

0.4 

- 0. 9 

6 

Double-Acting Bellows 

0. 02 

- 0. 06 

0.4 

- 0. 6 

6 

Aeronutronic BOA 

0. 13 

- 0. 17 

0.85 

- 0. 95 

6 

Proposed Lunar 

0. 15 

- Q. 18 

0. 02 - 

0. 04 

2 


Wheeled Vehicles 


By using the parameters given in the table, and assuming a vehicle 
density, y , nondimensionalized obstacle clearance and footprint pressure 
can be converted to absolute values as a function of vehicle size. Using the 
nomenclature h, b, and L for vehicle height, width, and length, this can be 
done as follows. The height is assumed equal to the width. 


M = y b" 


(r) 


S = Lb 


(i?b) 


(b) ( L S b) 
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This result is shown in Figure 79 for an assumed density of 10 lbm/ft . 
Typical vehicle densities are shown by Bekker5 , and 10 lbm/ft 3 is typical of 
small lightweight vehicles. Larger vehicles, especially typical military 
vehicles tend to be heavier than this. The result on Figure 79 can be ratioed 
by y 2/3 to give results for other densities. 

The obstacle clearance capability, h D , is given by 

h o = (r) fc) b 

■ (jf (it) (fJ (77) 

This is shown in Figure 80. 

An interesting presentation of power requirements of mobility systems, 
originally pointed out by Gabrielli and Von Karman, is derived from the work 
of Bekker and shown in Figure 81. The figure shows equivalent drag -to - 
weight ratio (sometimes called specific power) as a function of speed. The 
Gabrielli - Von Karman line, defined by D e / W = v /1720 in ft-lb-sec units, 
has been considered to be a limiting condition imposed by technology. 

Figure 82 shows the worm type of vehicles plotted to the same scale. It is 
seen that the estimated performance of the extension-contraction and rib- 
walking vehicles is in the area traditionally occupied by animals. The 
traveling wave is generally in the region of present cross-country vehicles. 
Although it is not anticipated that the traveling wave can match conventional 
vehicles (wheels) on a hard surface, it is anticipated that the range of soils 
in which it is useful will be much greater. 


MISSION SUITABILITY 


Mission Requirements 

The proper performance measures for a vehicle are strongly dependent 
upon the detailed requirements of the mission. Cost, the conventional all- 
encompassing measure of a transportation system, does not appear as a 
primary parameter in the evaluation of lunar vehicles. Parameters which 
might be considered are weight, size, reliability, power requirements and 
speed. Speed is not considered a primary requirement in the present spec- 
trum of missions, but will most certainly become more important in the 
future. Reliability is often the most important requirement of all. This may 
take at least two forms in the type of vehicles being considered here: reli- 

ability as a shelter, and reliability as a transportation system. This distinc- 
tion allows a tentative classification of the concepts into mission categories. 
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* Power requirements are generally an important consideration only in that 
they usually directly affect the size and weight because of the requirement for 
batteries or other energy storage or collection devices. Size and weight while 
operating are unimportant, but during the trip to the moon are very important. 
On larger vehicles some kind of folded, collapsed, or disassembled mode is 
required to meet the size constraint. 

The vehicle performance parameters must specifically consider the 
ability of the device to perform satisfactorily in the lunar mission environ- 
ment. Conditions that potentially pose design and operational constraints on 
a vehicle are high vacuum, temperature extremes, solar radiation, abrasive 
surface material, and meteorite bombardment. 


Ability to Fulfill Mission requirements 

The design studies in this report demonstrate some of the possible 
mission applications of these mobility concepts. The mobile shelters may be 
thought of first as a shelter, and must have high reliability of functioning as 
such. The ability to move, on the other hand, may be useful, but not vital 
in the performance of the mission. This reasoning tends to favor the 
extension- contraction device, which requires a minimum of mechanisms, 
moving parts, and appendages to the basic shelter structure to provide 
mobility. The reliability of this method of propulsion is no better than the 
confidence that the soil characteristics are as predicted and are fairly uni- 
form and reasonably smooth. Generally, this device cannot back out of a 
corner, nor can it move at all in a soil for which is was not designed. The 
extension- contraction vehicle requires special consideration in this respect, 
since it has been proposed that solar energy might be used directly for pro- 
pulsion, thereby avoiding the usual weight penalty for power consumption. 
Conversely, the lunar supply vehicle mission may require high reliability 
as a transportation device to avoid leaving the occupants stranded in a remote 
area, while its function as a shelter may or may not be important. This 
concept favors the traveling wave vehicle which claims, as its primary 
attribute, the ability to move on an extremely wide range of soil types. It 
also has the ability to bridge crevices having widths approaching one-third 
its length. The unmanned missions generally require mobility somewhere 
in between these two, and are possible applications of the double-acting 
bellows and rib-walking concepts. 

All bellows concepts studied can have all moving parts internally located 
in a protected pressurized environment to protect rubbing or sliding surfaces 
from cold welding in the extreme vacuum environment of the moon. 

All concepts are suitable for delivery in an undeployed length generally 
less than one- third their nominal length. A possibility exists of using bel- 
lows contraction as an emergency measure to provide additional protection 
from radiation during solar flares. If the bellows were made of metal and 
were sufficiently contracted, the effective thickness for shielding would be 
substantially increased. The degree of such protection which could be ob- 
tained was not established in the present study but some observations can be 
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FIGURE 81. EQUIVALENT DRAG-TO-WEIGHT RATIO 
OF PRESENT SYSTEMS 
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FIGURE 82. EQUIVALENT DRAG- TO- WEIGHT RATIO 
OF WORM-TYPE VEHICLES 
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made regarding the degree of contraction possible with metal bellows. If a 
welded bellows cross section is employed, such as shown in Figures 54 and 
55 for the lunar shelter, the bellows can contract nearly to their flat plate 
thickness and will be limited by other structural elements such as cleats or 
internal linkage systems. In the case of the large shelter it was estimated 
the shelter could be retracted to nearly one- sixth its maximum length. If 
a formed corrugated metal bellows cross-section is employed, however, 
which has sizable radii instead of cusps at the inner and outer edges of the 
convolutions, the length can be changed only about ten percent without ex- 
ceeding the fatigue limit of the material. 


Versatility Versus Complexity 

Of the mobility modes discussed here, the traveling wave mode is 
potentially the most versatile and efficient but probably will require the most 
complex actuating mechanism. The extension- contraction mode is the least 
versatile and efficient, but the simplest. Its low efficiency may be partially 
compensated for by its potential ability to directly convert solar energy to 
power for propulsion, although the practicability of this concept has yet to 
be established. The versatility of the double-acting bellows version of the 
rib-walking mode lies somewhere between the traveling wave mode and the 
extension- contraction mode. While it will require a relatively simple 
mechanism, it is limited in obstacle climbing ability and by roughness of 
ride. 
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CONCLUDING SECTION 


CONCLUDING REMARKS 


An analytical study has been carried out concerning the feasibility of the 
bellows vehicle concept for use on the moon. Since most of the work is based 
on motions which have been successfully used by animals, many of the re- 
sults have application to mobility in general and are therefore of interest in 
connection with other programs. 


Of the several modes of bellows motion studied the traveling wave mode 
appears to be potentially the most versatile since, with the proper phasing 
of adiacent bellows segments, it is capable of moderate to high speeds over 
rough terrain with a smooth ride and good maneuverability. Compared with 
current cross-country vehicles it has good propulsive efficiency, better 
obstacle climbing ability, and better ability to bridge crevices. It can be 
designed with a very low footprint loading suitable for use over a wide 
variety of surfaces including spiny rock, low bearing strength materials, 
and even fluid-like soils through which it is capable of swimming partially 
or completely submerged. It is adaptable for a pressurized membrane-type 
construction which not only reduces weight but also allows the vehicle lower 
surface to mold itself to the terrain and thus minimize vehicle bounce and 
lift-off under low gravity conditions. Its feasibility is dependent upon the 
availability of a suitable membrane cover material which will remain flexi- 
ble and abrasion resistant in the lunar environment. 


The rib-walking bellows motion where adjacent ribs contacting the ground 
are 180 degrees out of phase requires a large path radius motion to yield ac- 
ceptable drag-to-weight ratios. A double-acting bellows concept was devel- 
oped to provide such motion with only a moderately rough ride. This concept, 
when used in a brief design study of a small unmanned vehicle of modest range 
appeared quite suitable with regard to simplicity, ease of deployment, weight, 
and power requirements but was somewhat limited in ability to clear 
obstacles. In this study, feet were utilized to provide roll stability and to 
prevent contact of the metal bellows surface with the ground, thus keeping 
the structural design within the present state-of-the-art of bellows technology. 
Other possible configurations include side-by-side bellows and double-lobed 
membrane bellows. 


The extension- contraction bellows motion would appear to have limited 
use because of the high drag-to-weight ratio resulting from the sliding action 
of the lower surface over the terrain even with the use of a low-friction type 
coating and a variable geometry profile such as flexible cleats which retract 
during the forward stroke. The extension- contraction bellows configuration 
is, however, uniquely suited for incorporation of a novel solar -mechanical 
energy transfer system which from exploratory calculations holds promise 
as a lightweight, inexpensive means for supplying the mechanical energy 
for propulsion. A brief design study of using the extension- contraction 
bellows concept to propel a lunar shelter restricted to forward motion at 
very low speed over a short distance did not reveal any unusual problems 
regardless of shelter size. 
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All bellows concepts studied can have all moving parts internally located 
in a protected pressurized environment. All concepts are suitable for delivery 
in an undeployed length generally less than one-third their nominal length. 

In carrying out the present study fundamental analytical treatment was 
given: to a stepped surface planing over soil; to wave motions which might be 

used by a variety of vehicles and which included traveling waves in soft soil 
with no hysteresis; and to membrane bellows structural properties and shear 
transfer. It was also necessary to establish a lunar soil model more realistic 
than the ELMS model for studies involving short stroke and cyclic motion. 


RECOMMENDATIONS FOR FUTURE WORK 


For all the vehicle concepts discussed in this report, the principal problem 
area is the interface between vehicle and soil. The predicted performance of 
the various propulsion mechanisms should be checked by carefully planned 
tests supported by additional analysis and evaluation. 

The traveling -wave configuration which has been presented requires the 
development of an abrasion-resistent elastomer which will not deteriorate in 
the lunar environment. Further study should be undertaken to determine 
whether such a material can be produced. 

The performance of the various systems over rough ground has not been 
investigated. This performance, particularly in the case of the traveling wave, 
is strongly dependent on the degree of flexibility and control built into the body. 
Dynamic studies and tests are needed to obtain quantitive design information 
on these relationships. 

In order to obtain an overall evaluation of these concepts, mission-oriented 
trade-off studies should be performed, and resulting vehicle designs compared 
with other proposed vehicle designs. 
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UNITS AND CONVERSION FACTORS 


In most cases, the text and equations of this report use standard English 
units. The inch is used widely as the unit of length since the ELMS model 
was given in this form, and since inches are traditionally used in structural 
work. The only unusual unit used is pounds mass as the unit of mass. 

A pound mass is defined as the mass which weighs one pound on earth. In 
this system, a body's mass and weight are numerically equal on earth. 

The conversion factors to convert English units to the International 
System of Units are given below. The International System is abbreviated 
SI for Systeme International and is based upon the meter, kilogram, and 
second. These units are described and the conversion factors quoted here 
are given in Reference 7 and elsewhere. The first two digits of each con- 
version factor represent a power of 10 multiplier. For example, in the first 
conversion below, +03 1. 06 means multiply by 1. 06 x 10^. 


TO CONVERT 

TO 

MULTIPLY BY 

British thermal unit 

ioule 

+ 03 

1. 06 

foot 

meter 

- 01 

3. 048 

horsepower 

watt 

+ 02 

7. 457 

inch 

meter 

- 02 

2. 54 

pound (pound force) 

newton 

+ 00 

4. 448 

pound mass (lbm) 

kilogram 

- 01 

4. 536 

mile (U. S. statute) 

meter 

+ 03 

1. 609 

ton (short, 2000 lb) 

kilogram 

+ 02 

9. 072 

ACCELERATION 

foot/ second^ 

meter/ second^ 

- 01 

3. 048 

free fall, standard 

meter / second^ 

+ 00 

9. 807 

AREA 

foot^ 

meter^ 

- 02 

9. 290 

inch^ 

. 2 
meter 

- 04 

6. 452 

DENSITY 

3 

Ibm/inch 

3 

kilo gram /meter 

+ 04 

2. 768 

lbm/foot^ 

3 

kilogram /meter 

+ 01 

1. 602 
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TO CONVERT 

TO 

MULTIPLY BY 

ENERGY 

British thermal unit 

joule 

+ 03 

1. 06 

foot pound 

joule 

+ 00 

1. 356 

POWER 

foot pound/second 

watt 

+ 00 

1. 356 

horsepower 

watt 

+ 02 

7. 457 

PRESSURE 

pound /foot^ 

newton /meter^ 

+ 01 

4. 788 

pound /inch 

newton /mete r^ 

+ 03 

6. 895 

SPEED 

foot/second 

meter /second 

- 01 

3. 048 

mile /hour (U. S. statute) 

meter /second 

- 01 

4. 470 

VOLUME 

foot 3 

meter 3 

- 02 

2. 832 

• v.3 

inch 

meter 3 

- 05 

1. 639 
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APPENDIX B 


THE ENGINEERING LUNAR MODEL SURFACE (ELMS) 

The ELMS document, ^ is an attempt to present a uniform model of the 
lunar surface, in order to make the mobility studies of various contractors 
consistent with one another. It defines the mathematical form of the descrip- 
tion of soil properties, and specifies the range of soil property parameters 
that should be considered, as well as giving typical lunar surface profiles.. 

The soil description portion is based largely upon the work of Bekker, ^ 
and, since it is pertinent to the present study, a summary is given here. 

The relationship between footprint pressure, O’, and soil penetration 
depth, z, is given by 

17 ■ fr + K „) 

where b is the width of the print, and K c , K^, and n are soil parameters 
whose significance is described elsewhere in the List of Symbols. For the 
lunar surface, K c is taken as zero, and the values for K<p and n are given 
in Table B.l, together with a corresponding description of the soil and the 
expected description on the lunar surface. 

TABLE B. 1 


LUNAR SURFACE PENETRATION 


Slope 

Percent of Traverse 

Soil Constants 

Soil Description 

Angle 

Maria 

Uplands 

n 

K <p 





Lb 

(No 

Dimensions) 


Degrees 

Percent 

Percent 

iBH 





1 . 0 


Extreme soil 

0-4 

84. 0 

57.1 

0. 5 

BBH 

Loose, fluffy, dust 

5-7. 5 

10. 5 

23. 4 

0. 75 

1.0 


10-20 

5.0 

18. 3 

1.0 

3. 0 

Loose, dry, sand 

22-25 

0. 3 

1. 0 

1. 25 

6. 0 

Hard, compacted, sand 


0. 2 

0. 2 

- 


Hard surface 
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The "extreme soil" has been predicted by some theorists, but no specu- 
lation as to its distribution is made. The result of Equation (B . 1) with the 
range of parameters given in Table B. 1 is plotted in Figure 2 of this report. 

The shear stress, T, is given in its general form by 




T = (c + a tan <p) f(j) 

(B. 2) 



f(j) = Y 




y max 

(B . 3) 

y = exp I 

(-k 2 + V k 2 - 

V K ljJ - ex p[(- K 2 "V K 2 - 1 ) K lJ J 

(B.4) 


The function f(j) gives the variation of tractive force with j, the displace- 
ment of the soil from its initial location, and has a maximum value of 1.0. 

The shape of this curve is determined by soil parameters Ki and K The 
level of the maximum friction coefficient is given by the coefficient of cohe- 
sion, c , and the angle of friction of the soil grains, <p. The lunar surface is 
assumed to be a noncohesive material and c is therefore taken as zero. The 
friction angle, <p, is assumed to be 32 degrees which corresponds approxi- 
mately to dry sand, except for the "extreme soil" which is assumed to have a 
coefficient of <p =20 degrees. Slippage factors, Kj and > are conserva- 
tively chosen, low values reflecting a pessimistic prediction of vehicle per- 
formance. K\ generally ranges from 0.3 (average soil) to 1.0 (coherent 
soil), and from 1.0 (coherent soil) to 3.0 (loose sand). NASA analyses 
are based on Kq =0.2 and = 1. 25, thus yielding a conservative vehicle 
performance estimate. 

Using these values, Equation (A. 2) becomes 

r - tan <p / -0. lj -0.4j j 

- V / (B- 5 ) 

This relationship is plotted in Figure 1 of the report. 
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APPENDIX C 


MOBILITY OF EXTENSION-CONTRACTION TYPE OF VEHICLE 


As discussed in the body of the report, three soil models are considered 
as follows: 


Case I - — constant 

o 


T 

— (aft direction) = k 
Or a 


7* 

(forward direction) = 


Nominal values : k = l.Z L = 0.75 

a I 


k f = 0.625 


<T 

Case II - -jp proportional to displacement 


-j- (aft direction) = k^j 


7 

(forward direction) = k^j 

Nominal values : k = 0.475 in. 

a 


-1 


k f = 0. 396 in. _1 


Case III - ELMS Model 


4 = 4 ' 


O.lj _ e -0 . 4j 


0 


Nominal values : Aft C =1.58 

a 

Forward = 1.32 


The nominal constants given make the T /(7 for Case I equal to the peak value 
of the ELMS model, and the slope for Case II equal to the slope at the origin 
of the ELMS model. 
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The vehicle length is taken at any time as L (1 + e), where L is the 
minimum length during an operating cycle, and e is the vehicle strain. The 
length of the stroke is the Le. Since the constants given for the ELMS model 
are given for inches as the unit of length, inches will be used throughout this 
section. The motion of the vehicle will be referred to the midpoint, and since 
the vehicle is considered to have a uniform density, this will also be the 
center of gravity (eg). The point on the vehicle which is instantaneously at 
rest with respect to the ground will be labeled p. The vehicle weight, the 
tractive force, and the portion of the minimum length, L, on either side of p 
will be subscripted a and f for aft and forward of this point. At this time no 
external or inertial forces are considered. 


Analysis of Case I 


The geometry for this case is shown in Figure C-l. 
relationships apply 

The following 

H = kas 

(C. 1) 

W 

a " L (1 + e) 

(C. Z) 

L 

W = W ~ 
a L 

(C. 3) 

ii 

£ 

(C. 4) 

£ 

4- 

cti 

£ 

II 

£ 

(C.5) 

L = L + L, 

t / \ 



T 

a 



FIGURE C-I. CASE I T/a CONSTANT 


Since no external forces are considered, the location of the stationary 
„ point, p, is given by setting the thrust equal to drag (H a = Hf), and by sub- 
stitution of Equations (C. 1) through (C. 6) gives 


L a = k f L f 

(C. 7) 

k f 


L a = L IT-+T7 

(C. 8) 

a f 

k 


L f = L k + \ 

(C. 9) 

a 1 


The point p remains fixed throughout one extension or contraction stroke. 
The location of the eg with respect to p is 


d h = + e) 

a f 


and the eg travel during a half cycle is 


k - k- 

s = d, - h, = 1/2 Le , 

h ho k + k 


f 


(C. 10) 


<C. 11) 


137 



The work required to move this distance may be found as follows (where 
only the aft end of the vehicle is considered): 


dx 


dH a = k a W a L (1 + e) 


on element dx 


ds = 


x 


L ( 1 + e) 
a 


d(Le) for an incremental expansion d(Le) 


Work) 


-L e X (1 + e) 

J dH a ds 

o o 


t 


= 1/2 k W L e 
a a a 


(C. 12) 


Writing a similar expression for the forward end of the vehicle, and substi- 
tuting from Equations (C. 3) to (C. 6) gives the work for the whole vehicle 


Work) (ot = 1/2 WLe 


(C. 13) 


f 


Then, using Equations (C. 11) and (C. 13), the equivalent drag-to-weight 
ratio for the vehicle is 


D 

e 

“w" 


Work) 


tot 


Ws 


k k r 
a f 

k - k 
a f 


(C. 14) 


Analysis of Case II 

Case II is illustrated in Figure C-2. Equations (C. 2) through (C. 6) apply 
to this case, but the tractive force is given by 


H = 


t* 


(1 + e) 


T dS 
a 


W 


bL (1 + e) 0 


r L (1 + e) 

/ 


b k 


a 1 + e 


dx 


1/2 W k L e -f. 

a a L 


(C. 15) 
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(C. 16) 


H 


f 


1/2 W k f 




FIGURE C-2. CASE II T/cr PROPORTIONAL TO DISPLACEMENT 
Proceeding as for Case II, the following relationships result; 

H = H r 


k L 2 = k, L 2 
a a f f 




h = 


^F\ + V^7 
V*T 

Vk + V^7 


(C. 17) 


(C. 18) 


(C. 19) 


s 


1/2 Le 


ViT + V*7 


(C. 20) 
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dH 


- W k e 

a bL ( 1 + e) al + e 

a 


x b dx 


ds = 


d(L e) 


L (1 + e) ' a 

a 


Work) = 1/6 W k 

a a a a 


Work) tQt = 1/6 W(Le)‘ 


k r k 
f a 


(v^ + w) ; 


D 


k k r 

e i # *2 T a i 

W = 1 Le k - k. 


(C. 21) 

(C. 22) 
(C. 23) 


Analysis of Case III (ELMS Model) 

For this case, the point p is at rest only momentarily, and its location 
changes continuously. The variation of the location of this point has not been 
solved. However, it is known that at the beginning of a stroke the location of 
p is given by Case II, and it has been shown graphically that its departure 
from this location is not great up to a stroke length of about ten inches. 
Therefore, an approximate solution may be obtained by assuming the p 
remains fixed as given by Case II, and computing the resulting power 
consumption. 

The forces on a vehicle segment with one end fixed are shown in 
Figure C-3. The tractive force on the vehicle may be found as follows: 



Le 

/ U) d J 

0 


H = W 


T_ 

O ave 


(C. 24) 


H = 



( 


7.5 - 10 e 


■ s f 4-- u -■••*)* 

) (C. 25) 


0 

• 0.1 Le 


+ 2. 5 e 


-0.4 Le 
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FIGURE C-3. CASE III ELMS MODEL - VEHICLE FIXED AT ONE END 


The work consumed during a half cycle is given by 

,(1 + e) 


Work 


■ft 


dH da 


where 


, W T 

dH = bL (1 + e) "^ X) b dx 


da = 


d(Le) 


(C. 26) 


L (1 + e) 

Making the substitutions, and performing the first integration gives 

Work = W Fj(Le) 

where the function Fj(Le) is given by 


(C. 27) 
(C. 28) 
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Fj(Le) = 


r 


+ 2. 5‘ 


c 

r.io 2 ( + 

(Le) 2 

L \ io 

/ Le 

,\ -0. 4 Le 

\Ts + 

1 e 


-0. 1 Le 


+ 93 


• 75 ] 


d(Le) 


(C. 30) 

This function may be evaluated numerically, and the result is shown in 
Figure 9. 

For a complete vehicle which is free to move, the location of p is 
assumed to be the same as given for Case II, and the distance traveled per 
stroke is given by Equation (C. 20). The work is found by applying 
Equation (C. 29) to each end of the vehicle. 

Work) tot - TT WF l (L a e > + L WF 1 < L f e > (C. 31) 


where the function, F i is defined by Equation (C. 30). Then, by using 
Equations (C. 20) and (C.31), the equivalent drag-to-weight ratio is given by 


D 

e 

W 


Work 

Ws 


_2_ Vfl 

Le yc; - 



F 1 <V> 


+ Fj (L f e) J 

(C. 32) 


where L a and Lf are given by Equations (C. 18) and (C. 19). It can be shown 
that performing a first order approximation for small Le to this expression 
gives the same result as Equation (C. 23), the solution for Case II. 


Results 


The results of the analysis are shown and discussed in the body of the 
report. The work required to extend a vehicle which is fixed at one end is 
shown in Figure 9 , and is obtained from Equations (C. 12), (C. 21), and 
(C. 29). This much of the solution is exact within the limitations of the 
assumptions made. The equivalent drag-to-weight ratio is given by 
Equations (C. 14), (C. 23), and (C. 32), and is shown in Figure 10. Here, the 
additional assumption is made that the point p remains fixed at the location 
given by Case II, an assumption which is not valid for large stroke lengths. 
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APPENDIX D 




MATHEMATICAL MODEL OF WAVE MOTIONS 


The following is a mathematical description of motion which might be 
used by various types of vehicles. It is hoped that this model might prove 
useful in future studies of the mobility and dynamics of vehicles using these 
principles. The general type of motion which is to be considered consists of 
traveling waves along the lower surface of an elongated body, as shown in 
Figure D-l. 



If it is assumed that the vehicle length L contains N similar waves, the 
fundamental wavelength will be 


X - L/N w (D. 1) 

The vertical and horizontal components of the motion will be taken as sinu- 
soidal. This does not mean, however, that the wave shape is sinusoidal. In 
fact, the typical traveling wave vehicle mechanized by a series of cranks as 
discussed in this report is conspicuously nonsinusoidal in shape and actually 
is a curtate cycloid. 
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The coordinate system will be chosen as shown in Figure D-2 (fixed in 
vehicle). The ground plane will be represented by Z =0. 



FIGURE D-2. COORDINATE SYSTEM 


All points will be assumed to move in planes parallel to the XZ plane. If X is 
the distance in feet from the nose, and v w is the wave velocity in the aftdirec- 
tion in feet per second, the components of displacement of a point (x, y, z) 
with respect to the point (X, Y, Z) will be: 

Longitudinal Displacement: 

£ = A cos ~t- — (X - v t) (D . Z ) 

x x A ' w 

Transverse (vertical) Displacement: 

£ - A cos (X - v t + X ) (D. 3) 

z z A w o' 

In these equations, A and A are shown as constants, although theymight 
vary with Y or Z for certain more complicated types of motion. In order for 
useful thrust to be produced, Xo will normally be close to A/4, thus producing 
circular or elliptical local displacements. Values 0< X Q < 17 give propulsion 
in the forward direction. If X D is outside this range, propulsion in the oppo- 
site direction results. The displacements are shown in Figure D-3. 
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FIGURE D-3. DISPLACEMENTS 


One additional parameter is the number of frames or bulkheads per wave- 
length. This will be called n^. 

From Equation (D. 2), the longitudinal displacement velocity of a point x 
is 


i. 


x = 




2 IT v. 


w 


• 2 ir 
sxn (X 


v t) 
w 


which has a maximum value 


(D.4) 


v t 


m 


= 21TA 


w 


X X 


2 V f A 


This is the velocity with which the vehicle would move under ideal condi- 
tions on hard ground with no slip. Under these conditions, the speed of the 
wave relative to the ground (in the aft direction) is: 


v = v t 


m 


( 21TA * \ 
" V= y ~ X l ) 


(D. 5) 


That is, for A x = A/21T, the wave is stationary with respect to the ground 
(the centipede case). 
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Special Cases : 

The extension-contraction motion discussed in this report is a standing 
wave obtained from Equations (D.2) and (D. 3) by setting A = 0 a £ z and 


Putting 


A 

x 

“ 2 


cos - v w t) cos^-(X+v w t) 


( f Zir 

MS T 


x 


cos 2 TT f t 


■) 


X = -j- + 6, where 6 < <A gives: 

t a • 2 IT 6 _ _ , . 

£ = - A sin — r — • cos 2 TT ft 
x x A 




(D.6) 


2tr A 


• cos 2 7T ft. 


(D.7) 


(where f = frequency in cycles/sec) 

The rib -walking motion is approximated by setting 

X Q = A/4; n f = 2 

This gives for successive ribs, X = r* "y" (r = 0, 1, 2, etc.) 


£ = A cos 

* X X 




= A cos (r IT - 2 IT f t) 
x 


or 


= - A cos 2 IT £ t 
x 


+ A^ cos 2 TTf t 


(r odd) 

(r even) (D. 8) 


and 

£ -A sin (r TT - 2 TT f t) 
z z 

= - A sin 21T ft (r odd) 

z 

or = + A sin 2 77 f t (r even) (D. 9) 

z 

That is, successive ribs describe elipses, 180 degrees out of phase with 
each other. 



The traveling wave motion may be represented by setting 
X = 2 L, A x = A z = A, X Q = A/4 


giving : 


£ = A cos - 7 - (X - v t ) 

x L w ' 

£ = A sin 7 — (X - v t) 

z L w 


(D. 10) 
(D. 11) 


APPENDIX E 


MOBILITY OF A TRAVELING 


WAVE IN SOFT SOIL WITH NO HYSTERESIS 


Introduction. 

An analysis of the mobility of a traveling wave in soft soil with one set of 
assumptions regarding geometry and soil behavior has been completed. A 
wave made up of rotating cranks as shown in Figure E-l is considered. How- 
ever, the crank radius is assumed small compared to the wave length, re- 
sulting in a sinusoidal shape and certain other simplifications in the analysis. 
Many cranks are assumed, giving approximately a continuous wave. The soil 
normal force is assumed to depend only upon the depth below the original sur- 
face of the soil, and to be directly proportional to that depth. The soil fric- 
tion is assumed to be directly proportional to the normal force, i.e. , the 
classical friction coefficient. 

The principal results of the analysis are shown in Figures 22, 24 and 25 
in the report. Figure 22 shows the ratio of the vehicle forward velocity, V, 
to the crank velocity, Aoa Figure 25 shows the depth that a wave will pene- 
trate into the soil as a function of the planform loading. Figure 24 gives the 
equivalent drag -to -weight ratio, which is actually a measure of vehicle power 
consumption under the specified conditions. 

Analysis - Entire Wave Below Soil Surface 

Assume a traveling wave on the lower surface of a vehicle, which may be 
made up of many cranks as shown in Figure E-l. Each crank has radius A, 
and is rotating counter-clockwise with angular velocity co, and the wave is 
moving from left to right. 


See footnote on page 1. 
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FIGURE E-l. TRAVELING WAVE MADE UP OF CRANKS 
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If there are many cranks, and the radius is small compared to the wave 
length, then this will approximately represent a smooth, continuous, sinu- 
soidal curve as shown in Figure E-2(a). 



(a) SURFACE SHAPE 



(b) HORIZONTAL VELOCITY 

F05473 U 

FIGURE E-2. CONTINUOUS WAVE SURFACE 


The coordinate system of the traveling wave is fixed in the vehicle. For 
this case the entire wave will be assumed below the ground surface. 

The soil is assumed to have a linear relationship between sinkage depth 
and footprint pressure 


a = K z (E.l) 

<P 

and this relationship is assumed to hold even though the surface of the wave 
may be receding from the soil. The depth of the wave in the soil as shown in 
Figure E-2(a) is: 

z = Zq - A sin 277 ~ (E. 2 ) 
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The horizontal velocity of each point of the wave is: 


x = - Aco sin 2 IT - (E.2a) 

as in Figure E-2(b), where positive velocities are to the right. With respect 
to the vehicle, the ground is moving to the right with velocity V. Therefore, 
in Figure E-2(a), the portions of the surface between C and B are moving to 
the right faster than the ground, and are providing propulsion. All other 
points are contributing drag. 

The footprint pressure at any point is, from Equations (E.l) and(E.2): 

a = K (p Z 0 -KpAsinZTT f- (E.3) 

A constant friction coefficient is assumed so that 

T = jU a (E. 4) 

which is independent of velocity. Therefore, the propulsive thrust per unit 
of vehicle width for one wave length 

H = Jr d S 
p *'L 


= M 


K z 


cp o 




(E. 4a) 


The horizontal component of the force parallel to the surface is taken here by 
using the area of the surface which is projected on the horizontal plane; i.e. , 
dS = dx. Similarly the drag force on the wave surface is given by: 




X 

T 


T d£ 


x 


TdS 



T dS 


0 


2 


f X TdS 
J B 


(E.5) 



The net propulsive thrust on the wave is H - H^, and from Equations (E. 4a) 
and (E. 5), 


H 


net 


B 

/ 


TdS 


/ 


TdS 


\ / -o ~ \ 


-J Td s 

(E.6) 

B 


21T -J- - cos 2TT 

(E.7) 


However, from inspection of Figure E-2, it may be seen that: 

2(X - B) 


b -c ; 


-j B « C 

-cos 27T — = cos Z7T — 

and Equation (E. 7) may then be written as: 

— - _ 4 (^ _ jj\ + — — — • 4 cos Z7T — 

UK z U ' z„ 2ff X 

p <p o 0 


(E.8) 


If there are other forces acting on the vehicle which total D per unit of 
vehicle width for each wave length, then: 


H - H , = D 

P d 


r-> 


Equation (E. 8) can then be written as 

D 


0 

This may be written as: 


(iK + 4(X ' B) = 4 z„ 2V — X 


A X r B 

cos 2«- 


(E.9) 


(E. 9a) 


U , B' 1 A 


2TT z, 


cos 211 


(E.10) 


where the following terms are defined: 

D 


U = 


z o x 

B' = X - B 


(E. 11) 
(E. 12) 
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and where 


0 _ B ^ - B ! 

cos 277 = cos 2 77 


(E. 13) 


A graphical solution to Equation (E. 13) is given in Figure E-3. The 
intersection of the two families of curves on Figure E-3 is the solution of the 
equation for B'/X for that set of parameters. 


From Figure E-2(b), the ground velocity is: 

B ' 

v - sin - j- • 2 ^ 


(E. 14) 


0.16 


0.6 


0.12 


0.08 


0.04 


0.4 


U = 


0.2 


jU K,z \ 
0 o 


20 


40 


60 


360 (DEGREES) 

A 


- 0.2 


2 

* 

/ - 
/ Z 

/ o 

/ 


■1 

MB 


/ 

y-0.6 

|9 


i 

Wa 

/ 

,-0.8 


MB 

n 

fM 

H 


80 


100 


F05472 U 

FIGURE E-3. TRAVELING WAVE POINT OF ZERO VELOCITY 
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This result is plotted in Figure 22, by using the solution for B'/\ given by 
Figure E-3. In this figure, U is designated by the terminology which will be 
given by Equation (E. 17) and (E. 19) below. For clarity of presentation, the 
depth parameter is converted to z m»the maximum depth of the wave below the 
surface (rather than mean depth, Zq). This is given by: 


z 

m 

A 



(E. 15) 


Therefore, in Figure E-3, the curves for z m /A from 2.0 to infinity represent 
the present case of a fully submerged wave. The curves for z^/A less than 
2. 0 are for a wave partially above the surface; this case is discussed later. 


The vehicle weight (one wave length one unit wide) is, from Equation (E. 3), 


W, 


i 


K z _ - K A sin 2 IT r-1 dx 

(p 0 <p XI 




* 0 X 


(E. 16) 


This is the same weight that would result if the vehicle had a flat bottom, at 
depth zq below the ground surface. From Equation (E. 16), the term U given 
by Equation (E. 11) may be written as: 


1 D _ 1 External Drag 
JI ’ " ju * Weight 


(E. 17) 


If the drag, D, is the gravity force when the vehicle is on a slope of angle c p, 


D = Wj tan 9 


(E. 18) 


and Equation (E. 17) becomes: 

U = — tan 0 

M 


(on a slope) 


(E. 19) 


The power consumed by one wave length of the wave one unit wide is: 

Power = Force x Velocity (E.20) 

The power to drive the cranks as shown in Figure E-l is due to forces parallel 
to the wave and forces normal to the wave. We have assumed, however, that 
the pressure, 0, is acting normal to the surface, and is the same when 
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the wave is receding from the soil as when it is advancing toward the soil. 
Therefore, the power used in overcoming the normal force component is re- 
covered in another part of the cycle, and only the frictional component need 
be considered. 

From Equations (E. 3) and (E.4), the frictional force on an element of the 
surface in Figure E-Z(a) is: 

dP f = T dA 

= jU K A z n [1 - — sin 277 ) dS (E.21) 

0 \ Z 0 X / 

The required area increment is the actual area along the wave surface and is 
given by: 


dS 


dx 
cos 8 


(E. 22) 


where 8 is the angle shown in Figure E-2(a), and tan 8 is the slope of the 
wave at location x. The slope is given by 


tan 8 - — A sin 277 

dx 


x 

A 


Therefore , 


= A 


2 77 

X 


cos 277 -j— 


cos^ 8 - 


1 


1 + tan 8 


cos 8 


■H* 


cos 277 



(E. 23) 


(E . 24) 


The velocity required for Equation (E.20) is the tangential velocity of the sur- 
face relative to the vehicle, and from Equation (E.2a), 


_ -A CO 

t " COS 8 ~ COS 8 


V, = 


sin 2 77 -r— 


(E . 2 5) 


If the vehicle has no forward motion, the frictionforce defined by Equation (E. 2 1) 
is always opposing the direction of the motion given by Equation (E.25). The 
following discussion leading to Equation (E.29) is restricted to this case. The 
more general case with forward motion is treated as a correction to the zero- 
velocity case and is illustrated in Figure E-5. 
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The power may now be written as 


✓ 



v t T 



(E. 26) 


Note that T is always a positive number, but that v t is negative in the first 
interval, and positive in the second interval of Equation (E.26). The division 
Ox tire integral with the sign shown in this equation then makes the power 
positive, and additive over the whole surface. 


Adopting the nomenclature 



f(x)dx + 



X 


f(x)dx 


2 


A 


/ 


f(x)dx, 


(E. 27) 


the power may be written by using Equations (E.22), (E.25), (E.24), and 
(E. 26) as, 



Performing the integration and combining terms gives: 

p i- 4 ** I Vo A “^r [ 1 + t( 2 " f) 2 ] (e. 28 ) 

If the wave is relatively flat compared to the wave length, the second term of 
Equation (E.28) becomes small, and the power may be written as 


p r -r'' K <p z o x A " < E - 29 > 

The power given by Equation (E.29) is the power to move the wave surface 
when the vehicle is not moving, and power is being put into the wave internally 
at all points. Figure E-4 shows the frictional force on a wave which is on a 
vehicle with forward motion. The cranks as in Figure E-l are rotating 
counter-clockwise. The plane A-A represents the level on the surface at 
which the wave velocity is equal to the vehicle velocity. The arrows below the 
wave represent the direction of the frictional force on the wave due to relative 
velocity between the wave and the soil. The arrows above the wave represent 
the direction of the velocity of the wave with respect to the vehicle. Over 
most of the wave surface, the friction force is opposing the motion, and the 
friction is, therefore, extracting power from the vehicle. However, from 
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A/2 to C, and from B to X , the reverse is true and the frictional force is 
putting power into the vehicle mechanism if that mechanism is reversible. 

The power represented by Equation (E.29) represents vehicle power con- 
sumption only for B = X (B 1 = 0). For B' t 0, the power given by Equation (E.29) 
is too large by two times the power represented by the wave in the increments 
X /2 to C and B to X , or four times the power represented by the increment of 
wave from B to X , if the mechanism is reversible. The power that the 
vehicle can extract from the portion of the wave from B to X is, from 
Equation (E.26), 

X 

XP = J v T ds (E. 30) 

B 

and the true power consumption for a reversible mechanism is 

P = Pj - 4 c.P (E. 31) 

From the step preceding Equation (E.28), Equation (E. 30) becomes 
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FIGURE E-4. FRICTIONAL FORCES ON A WAVE 


Dropping terms containing (A/X) 2 gives the approximate soluti< 


M K <p z 0 Aa> = 2 


£ [* - cos 2 ”? 
- J sin 2 • 211 


Then, combining (E. 29) with (E. 31) and (E. 32) gives 

p * ¥^ k o z 0 X (A “» f z (f • 


where the function f^ is defined as 


, /B ! A \ 0 R ! 1 A /„ B' 1-9 _ B \ 

f 7 I t — , 1 = cos 2 7T — - — — 1 2 77- — - t sin 2 • 277 — ] 

2 \* Z Q I X 2 z 0 y x 2 x/ (E. 

Using the power given by Equation (E. 33), the weight given by 
Equation (E. 16), and the velocity given by Equation (E. 14), the equivalent 
drag-to-weight ratio may be found by: 


fU Jl 

X ’ 


(E. 35) 
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The function is defined by Equation (E. 34), and B'/\ is given in 
Figure 22 of the report. The result of Equation (E. 35) is shown on Figure 24 
in the report. This equation applies only to the fully submerged case which 
corresponds to Z m /A from 2.0 to infinity. The penetration depth, z m , is 
defined by Equation (E. 15). 

Part of Wave Above Soil Surface 

A case of a wave which is partially above the soil surface is illustrated 
in Figure E-5. 


GROUND 

SURFACE 
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FIGURE E-5. PARTIALLY SUBMERGED WAVE 


The location of the intersection of the wave with the ground surface is 
defined as E and is given by 


. E -*0 
sin 271 — = — 7— 




by: 


Proceeding as for Equation (E. 6) through (E. 10), the net thrust is given 


* R 

z~ B 


H 


net 


L “ k 4 


MK sin 2Tf — + z^Jdx -2 


) dx - 2 1 


^ K <p\ A sin 277 \ + z o] dx 


o> 


H 


net 


A 


■z?^ cos ^ 


B' , . . . B* , , o_e\ 

£ It — -6 IU a C* II — T l CO S c* tt — j 


fo /x 

A 12 


2 B' - 2 B' + 2 E 

Setting this equal to external drag, gives 
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(E. 37) 


The vehicle weight is given by 
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W. = K 

1 <p 
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A sin 27T + Zq | dx - 2 
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W. = K AX | — cos 27T 

1 <p )n 


Combining Equations (E. 37) and (E. 38) gives 
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(E. 40) 


For zq/A = 0, this becomes 

cos 27T x = i + ^ 

which gives a closed form solution for the half- subme rged case. 

Equation (E. 39) can be solved graphically in the same manner as was 
done for Equation (E. 10) in Figure E-3. If the equation is first multiplied 
through by zq/A, Figure E-3 can be used to obtain the solution, since the 
left-hand side of the resulting equation is represented by the cosine curve for 
A/zq = 1. 0 on Figure E-3. The right-hand side of the equation is a straight 
line for any combination of zq/A and D/jUWj and may easily be sketched on 
Figure E-3. The solution is the intersection of the two lines. The equation 
i s valid for zq / A from -1.0 to +1,0. The result in terms of penetration 
depth defined by Equation (E. 15) is given on Figure 22. The curves for d/A 
less than 2. 0 are from Equation (E. 39); curves for d/A greater than 2. 0 are 
from Equation (E. 10). 


The power for the special case of a wave half submerged (zq = 0, or 
d/A = 1.0) may be easily obtained. In this case, the portion of the wave from 
X /2 to C and from B to X is absorbing power from the soil, and the portion 
from C to B requires a power input from the vehicle as shown in Figure E-4. 


Then, using the form of Equation (E. 26), 
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(E. 41) 


For this case, zq = 0, and Equation (E. 41) can be multiplied through by 
Zq/A, and the terms containing a zq are dropped before making the evaluation. 

Using the relationships 
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and dropping terms containing 


(A/X) 2 


gives 


where. 



(E. 45) 
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By using the power from Equation (E. 45), the weight from Equation (E. 38), 
and the velocity from Equation (E. 14), the equivalent drag-to-weight ratio 
becomes, 


D -n mf 

e _ P It 

~W ~ VW = f* T 



This is shown on Figure 23 as the line for z m /A = 1.0. 


(E. 47) 


Results 

If the amplitude of the wave is small compared to the wave length, the 
significant geometry parameter is seen to be the ratio of maximum soil 
penetration to wave amplitude, z m /A (or alternately zq/A as given by 
Equation (E. 15)). The relationship of z m / A to the vehicle planform loading is 
given by Equations (E. 16) and (E. 38) and is plotted in Figure 25. The 
velocity ratio is given in Figure 22, and the equivalent drag-to-weight ratio 
in Figure 24. These are given as a function of a parameter which may be 
interpreted either as a slope, or as the ratio of some external drag force to 
the vehicle weight, divided by the friction coefficient of the soil. 

In Figure 22, the maximum slope capability is determined by the point at 
which the velocity ratio goes to zero. This is plotted as a function of z m /A 
on Figure 23. It may be seen that a vehicle that is less than half submerged 
(z m /A < 1.0) can climb a slope equal to the soil friction angle. 
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APPENDIX F 


CORRUGATED BELLOWS METHODS 


Design Charts 


To define the limitations of the welded metal type of bellows, and to aid 
in vehicle preliminary design, a number of graphs are shown which are based 
upon the equations in the report. Two materials were assumed. Unless 
otherwise noted, the material is stainless steel with a fatigue limit of 
25, 000 psi, Young's modulus of 28. 5 x 10° psi, and mass density of 
0.283 lbm/in. . In a few cases, results are shown for fiberglass with a 
20, 000 psi tatigue limit, Young's modulus of 3 x 10&, and density of 
0.065 lbm/in. The difference between these materials is most noticeable 
in its effect on allowable deflection, since the largest difference in the 
materials is in Young's modulus. The gravity factor g/g is taken as 1/6. 


Equation (47) in the report body is plotted on Figure F-l for stainless 
steel. The ordinate of the graph is f(d/b), but this parameter is not required 
except for computation and the scale is omitted. This graph may be read as 
in the following example. Enter with the desired d/b at point A, follow the 
dashed line to the desired working pressure at B, and to the desired width, 
b, at C. The useful weight to gross weight ratio is read at D. For a dis- 
cussion of the dashed crossing lines, see Equation (F. 7) below. 


Figure F-2 is the same thing for fiberglass for two pressures. The 
change in material affects the curve on the right only by the change in stress 
level; the curves on the left are affected only by the change in density. 


The vehicle footprint pressure is defined as the gross weight in the local 
gravity divided by the projected area at the bottom surface. This is given by 



From Equations (44), (45), and (39) of the report body, this becomes 
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where f (x) is defined in Equation (48), This is shown for stainless steel in 
Figures F-3 and F-4. The footprint pressure is low for small pressures 
because the bellows material becomes light and cannot support much payload. 
At high internal pressures, the payload becomes small because bending 
stress in the bottom surface is used in supporting the internal pressure. The 
effect of width is seen to be small. 

Assuming that the vehicle cross-sectional area is equivalent to a semi- 
circle of diameter, b, the force exerted on the ends by the internal pressure 
is 


R = $ b 2 p. (F. 3) 

The maximum bellows spring force is, from Equation (36) modified for a 
section with circumference equal to b, 


P 

s max 


2 


b dp i 


1 ,2 d 

b rr p. 

2 b 


(F.4) 


From the Mobility section, the maximum force required to overcome ground 
friction is (incremental extension) 


P r = l(u + (M + M )L (F. 5) 

f 4 ' *a 'r g ' u e ' e 

& e 

where (M u + M e )L e is the vehicle gross weight fully loaded. The require- 
ment that the internal pressure at least be able to overcome bellows spring 
force plus ground friction may be stated as 


P. > P 
i — s max 


+ P 


Substituting in the above relationships, this becomes 
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If this is made an equality and solved for pj, the result is the design 
pressure, p cr it» which will just overcome spring force and ground friction 
in the vehicle which is designed to that pressure as a structural criterion. 
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Substituting the footprint pressure from Equation (F. i) gives 


P crit 2 b ^ 


+ 

a P f' IT - 


d 

b 


(F. 7) 


1 where p is given by Equation (F. 2). Note that footprint pressure, p, is a 
| function of pressure, so that p cr ^ now appears on both sides of Equation (F. 7). 
| Ignoring the small effect of W in Figure F-3, p cr it a function only of d/b 
and the parameter L e /b (JLt a + Mf)> and may be solved graphically. 

The result (for stainless steel) is shown in Figure Fr5. This may be 
interpreted as the minimum design pressure for given parameters, d/b and 
Le/b (M a + Pf)f which will result in a vehicle which can extend itself to 
i maximum allowable length by internal pressure alone. This same result is 
also shown on Figure F-l as the dashed crossing lines. 

Equation (38) may be solved for the number of convolutions, N 

T, 


N = 


(i-V5) 


(F.8) 


2d 


From Equations (34) and (35): 


max 


F 2 t 
3 p. E 

r i 



(F.9) 


The deflection of the vehicle is given by N 6 and by using Equations (F. 8) 
and (F. 9): 


N6 
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1 

nr 




(F. 10) 


which is the maximum vehicle strain referred to the maximum extended length. 
Note that this is a function of only the material properties and the pressure. 
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MINIMUM DESIGN PRESSURE (PSI) 


A VEHICLE DESIGNED TO A 
LOWER PRESSURE THAN GIVEN 
CANNOT EXTEND ITSELF TO 



b 

F05480 U 


FIGURE F-5. MINIMUM ALLOWABLE DESIGN PRESSURE 
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MINIMUM DESIGN PRESSURE (PSI) 



Equation (F. 10) is plotted on Figure 30 for both stainless steel and 
fiberglass. For the large strains that result with fiberglass, the two dimen- 
sional Equations (28) through (31) no longer apply and the result is not valid. 
However, the large difference between steel and fiberglass in allowable 
deflection is evident and is valid within reasonable limits. 

The bellows spring force at this maximum deflection is: 


s max 


1 , 2 d 
2° b p i 


(F. 11) 


This is plotted on Figure F-6, and is seen to be independent of the material. 
Combining Equations (45) and (39) gives the vehicle empty mass: 


M = 3 y b 2 £ • 
e b 



(F. 12) 


This is plotted on Figure F-7 for stainless steel. This quantity may be 
critical in a design because of the requirement for transport to the moon. A 
minimum value of empty mass which defines the limit of what can be achieved 
by using the design principles assumed here. This is obtained by using d/b 
at each pressure which will result in the allowable stress being caused by 
beam loading of the bottom due to internal pressure only. This design can 
carry no payload; nor can it withstand the bending stress due to its own struc- 
tural weight. This relationship is the solution of Equation (43) for M u and y 
equal to zero, giving 



(F.13) 


This gives the value of d/b at which the curves on the right in Figure F-l 
intersect the abscissa, and is also plotted in Figure 31. This will be re- 
ferred to as the minimum weight design. 


Using the value of d/b given by Equation (F. 13) in Equation (F. 12) gives 
the emptv mass of the minimum weight vehicle. This is plotted in 
Figure 32. 
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FIGURE F-6. BELLOWS FORCE AT MAXIMUM DEFLECTION 
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If the vehicle is operated empty, the footprint pressure is the empty 
weight per unit length divided by the width. From Equation (F. 12) this gives 


= 3 ^ y 


■(#) 


(F. 14) 


Using the d/b for the minimum weight vehicle in this expression gives the 
empty footprint pressure of the minimum weight vehicle. This is shown on 
Figure 33. The thickness of the bellows material required for this mini- 
mum weight design may be found from Equation (34) and is shown in 
Figure 34. 

For a vehicle of unit length, the enclosed volume is equal to the cross- 
sectional area which is given by 


(F. 15) 


as for Equation (F. 10). Using the empty mass from Equation (45), and 
equivalent thickness from Equation (39), then gives the empty mass to 
enclosed volume ratio 


M ~ „ , 

e . 24 d 

V 7T / b 


lr\IE 

4 \l p. 


(F. 16) 


This is shown on Figure 35. For comparison, the ratio is also shown for 
a cylinder and a sphere designed to contain the same pressure. These latter 
structures are much more efficient than the bellows vehicles, since all 
stresses are taken in tension. 


The Mobility section discusses the concept of a foot planing on a granular 
type of soil and gives an expression for the maximum allowable footprint 
pressure if this planing is to take place. This may be written as 



where jL is the length of a planing foot in the direction of motion, y g is the 
soil mass density, and R and J are functions of the soil friction and the 
angle of attack of the foot. If it is assumed that each corrugation of the 
bellows has a foot, the maximum effective length of each foot is given as for 
Equation (38) as 


i 

max 




+ 



(F. 18) 


when the vehicle is at maximum length. Then from Equation (F. 17), the 
maximum footprint loading for planing is 


P 


al 




(F. 19) 


Setting this equal to the empty vehicle footprint pressure given by 
Equation (F . 14) defines the conditions for a vehicle design which will just 
plane. This gives 



RJ 


where y is the density of the bellows material and 7 is the density of the 
soil. Solving this for p^ gives 

•fk = JR ’ ' + ^( 3 ^ JR " ') " 1 ] <F - 20) 

y 

For large, this becomes approximately 

vk JR ' ') ,F ' 2I) 

This relationship is shown in Figure 36. The pressure p^ given is the 
maximum design pressure which may be used if the feet are to plane under 
the stated conditions. 
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Calculation of Bellows Parameters 
for Large Lunar Shelter 


The following calculations are based upon both the preceding equations 
and the vehicle characteristics given in the Design section of the report. 
Alternately, most of the values could be read from the design charts. 

In the nomenclature used for the bellows analysis, the assumed values 
are as follows : 

b = 240 inches 

~ = 0. 045 

b 

d = (240)(0. 045) = 10. 8 inches 

p. = 2 psi 

L e = (60)(1 2) = 720 inches 

F = 25, 000 psi 

y = 0. 283 lbm/in^ 

The material thickness is 



= 10 


.8-yf 


2 

25, 000 


10.8 

TTT 


0.0963 


The equivalent thickness at the maximum extended length is 



The vehicle empty mass is 


L M 
e e 


(‘ + r) b ^e L e 

(2. 57) (240)(0. 283)(0. 374)(720) 


= 47, 000 lbm 


i he useful loud is the weight which cun be Tended into the vehicle floor 
without exceeding the allowable stress when the floor is supported as a 
simple beam supported at the ends. This is given by 


L 

e 


M 


u 


(720)(240) 


1.73 x 10" 


»• 


(25, 000)( 10. 8)(0. 374) 
(240) 

|*(2. 34 - 2)(6) - 0. lO^j 


(1. 73)(10) 5 (2.04 - 0. 11) 
334, 000 lbm 



(0. 374)(0. 283) 


The maximum gross weight is 334,000 + 47,000 = 381,000 lbm and the 
useful load to gross weight ratio is 


334. 000 

381. 000 


0. 877 


At maximum length, the footprint pressure on the moon is 

„ _ 381,000 1 _ n . 

P min (720)(240) ' Z °* 368 psi 

The maximum deflection of the bellows from the neutral position is 


6 

7 " 


F 2 t 

^Pi E 


(25, OOP) 2 (0. 0963) 
(6)(2)(28: 5 x 10 6 ) 


= 0. 176 in. /disc 
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The number of discs may be solved from the given maximum extended length 
which is based on an assumed 1:4 slope of the discs 


2 N 


L 

e 





720 


(10. 8)(0. 25 4- 0. 0089) 


257 discs 


Then the maximum vehicle deflection is 

6 N = (0. 176)(257) - 45. 2 inches 

The vehicle neutral length (nonstressed position) 

L = 720 - 45 = 675 inches 
e 

and it operates with a stroke of 445 inches about this point. The minimum 
length is 

720 - 90 = 630 inches 


and the maximum footprint pressure is 

381, 000 

p (630)(240)(6) 


0. 42 psi 


The spring force developed by the bellows at the maximum deflection is 

P = i (2. 57) b d p. 
s 6 i 

= (2. 57 )(240)(1 0. 8)(2) 

= 2220 pounds 

t 

The cross sectional area of the vehicle is 


1 1Tb 1 2 _ 1 IT ( 240 ) 2 

2 4 2 4 


22, 600 in. 


The force exerted on the ends by the internal pressure is 

P. = (2)(22, 600) = 42, 000 pounds 




| At the maximum deflection this leaves Pj - P s = 40, 000 pounds of force 
to overcome friction while moving. 

The shelter volume at its nominal length is 

V = L S 
e c 


= (22, 6 00 ) (67 5 ) 


= 15. 3 x 10° in? 

and the empty mass to volume ratio is 

47 x 10 6 0 ^ 3 lbm 

-r - 3. 07 x 10 ^ 

15;3x 10 b in. 

j The closed length is 


L = 2Nt 
c 

= (257)(0. 0963) = 24. 8 inches 

The density of cargo which would fill the vehicle with its maximum useful load 
(loading density) is 


■■■ 334 . »„ 000 (1728) = 37.7 

(15.3x10) ft 
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ABSTRACT 


A feasibility study bas been carried out of various 
concepts of using a bellows to provide cross country 
mobility, particularly on the moon. Effort was directed 
toward analytical treatment of bellows mobility and 
of bellows structure including both corrugated and 
membrane bellows, and applying the results to several 
brief design studies of typical missions. All concepts 
studied can have relatively low footprint loading, can 
have all moving parts internally located in a protected, 
pressurized environment and can be considerably 
foreshortened during delivery. 

The concept wherein the bottom surface of the 
vehicle moves in the form of a traveling wave similar 
to the motion used by snakes and centipedes appears 
to be potentially the most versatile and is shown to 
have performance superior to conventional off-road 
vehicles. The double-acting bellows concept, which 
is a variation of rib-walking, has good performance 
suitable for unmanned vehicles but does have some 
restrictions. The extension-contraction bellows 
concept, wherein propulsion is derived from asym- 
metry built into the lower surface is limited in 
performance but is uniquely suited for incorporation 
of a simple solar-mechanical energy transfer system 
which appears promising. 

Alcn included in fundamental analvtical treatment 
of a stepped surface planing over soil, of wave motions 
which might be used by a variety of vehicles, and of 
the structure of membrane bellows. A lunar soil 
model suitable for short stroke and cyclic motion 
studies was developed. 
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